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ABSTRACT 


Administrators are often confronted with problems for which 
there exist several distinct measures of success. Such problems can 
often be expressed in terms of linear programming models with several 
linear "criterion" functions instead of single objective functions. 

A variety of techniques have been applied to multicriterion problems, 
but there exists a need for one which does not assume technical soph- 
istication on the part of the decision maker and which provides valid 


toy 
CHE 


solutions with minimum effort. "Efficient Maniffold Presentation, ' 
approach used here, is based on the concept that the ideal solution 

to a multiobjective problem must be a Pareto optimal solution. In many 
cases simply narrowing the set of candidate solutions to the set of all 
Pareto optimal solutions may enable the decision maker to find the com- 
promise being sougat. A finite method for finding an expression for 
the set of all Pareto optimal solutions to a Lincar program with 
multiple linear criteria is presented. Two processes are involved; 
first, the discovery of all Pareto optinal vertices of the feasible 
region, and secondly a grouping of these into sets each of which 
defines a convex polyhedrcn of Pareto optimal possibilities. Alter- 
Nabe verenons of thersecend prevess-are=suggested ior use- under varying 
circumstances. As an example of the applicability of the method, a 
general framework for modeling enrollment and staffing policy in an 
educational institution is presented, and a speci.ic model is develosed 


for the School of Healtn Services or tre Johns Hoszins universitr. 
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CHAPTER 1 


INTRODUCTION 


1-1 PROLOGUE: THE DEAN'S PROBLEM 

A college dean is pondering the related questions of how 
many students his school should admit and how many faculty he should 
have available (or perhaps how many he should terminate) in prepara- 
tion for an approaching academic year. He can predict readily 
enough the effects of his decisions for this upcoming year and per- 
haps, with somewhat more difficulty, for the year following that, 
but he also realizes that the potential impact of his decisions may 
linger long beyond. He therefore decides that he would like to 
address his problems in the context of several years, and after 
careful thought concludes that he can accurately assess the implica- 
tions of five consecutive years of enrollment and staffing decisions. 
More specifically, he is confident that he can project mumbers of 
students and numbers of faculty; from these he can derive a variety 
of revealing measurements, such as numbers of graduates, annual 
budget surplus, student-faculty ratios, differences between current 
and projected enrollments, and so on. Each such measurement is sig- 
nificant and none of them subsumes the others. Perhaps the Dean 
makes several such projections, each for a different set of enroli- 


ment and staffing assumptions. He measures each projection according 
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to those criteria he finds meaningful. It is quite likely that the 
criteria do not all favor the same projection. Nevertheless the Dean 
is able to choose a preferred projection from among the several 
candidates, one for which the strengths of some criteria overbalence 
the wealmesses of others. But is this the best possible projection? 
The Dean may well wonder whether he has overlooked an even more 
favoraole possibility. He may also wonder how he would recognize 

the best possible projection if he were to stumble upon it. 

Suppose that the Dean is able to enumerate the various cri- 
teria by which he measures alternative projections. One desirable 
preperty of this list is that it should be complete, in the sense 
that the information contained in the measurements according to these 
criteria is sufficient for him to be able to choose from among any 
alternatives so measured. Another desirable property is that all 
other things being equal, a projection with a lower "score" accord- 
ing to one criterion should not be preferable to one with a higher 
score. This should remain the case regardless of whether the scores 
are high or low, and regardless of how great or small the difference 
is between the scores. One must assume, also, that the Lean is able 
to specify the restrictions within wnich his projections must remain. 

The general practice has been, and often still is, for the 
Dean to base his decisions upon scratch pad calculations. The gener- 
al technique is frequently as follows. First he selects a set of 
appesilingedecisions, projects their impact, and scores the results 


according to the various criteria. Casting a critical eye over these 
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seores, he adjusts the decisions ima manner which he hopes will 

ameliorate some of the scores’ deficiencies, and repeats the pro- 
jection and scoring. This continues until the Dean concludes that 
he can do no better, or until he gets tired, or runs out of time. 

If the Dean's institution has some sort of plenning office, 
it may be possible for him to shift the burden of making the pro- 
jections and computing the scores to a computer. The formulae em- 
ployed to make and score the projections may constitute one of the 
many familiar “input-output” models. There are several ready-made 
computerized input-output models available, some of them extremely 

etailed. But regardiess of tha devail and precision of such a 
model, its focus is on projection rather than on decision. The Dean 
must nevertheless continue to hypothesize the decisions and judge 
the results. No input-output model can guarantee that the decision 
maker will be able to find a feasible alternative which suits him 
above all others, and none can guarantee that he will be able to 
recognize such a point should he arrive at one. 

The Dean’s planning office, if it is capable of more sophis- 
ticated institutional research, may suggest another approach in- 
tended to lift from the Dean the burden of decision as well as that 
of projection. A superficially straightforward technique consists 
of expressing the problem as a "mathematical program", and then 
solving this mathematical problem to discover the best decisions. 


Ih 


D 


mathematical programming formulation might be as follows: 





Maximize 0/1/1209 £, (x) ] 
Subject to g, (x) 2 
B(x) £ 0 


4 © © 


በ 
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where x is a vector with components X;,X2,++-,%, representing the 
enrollment and staffing decisions made. The inequalities gxi(x) ፥ ዕ 
express the limitations or "constraints" within which the institu- 
tion operates. Each f(x) is a function giving the score that de- 
cisions x achieve with respect to criterion 1, and the function H, 
the problem's "objective function", numerically expresses the Dean's 
satisfaction with any projection measuring f, (x), fo(x), +++, f}.(x) 
with respect to the various criteria. The functions f; (x) and g; (x) 
are relatively easily expressed: they are, after all, the stuff of 
"nich the many input-output models are made. The satisfaction 
function H[f4(x),fo(x),...,fy(x)], on the other hand, is exceedingly 
elusive, for it involves the comparison of incommensurables--the 
apples and oranges, so to speak, of the educator's trade. Ths Dean 
is able to directly compare two alternative projections and to pick 
the one which satisfies him most, but creating a function H which 
when confronted with two alternatives will unfailingly select the 
one which he would choose, is a very difficult task. Of course, the 
institutional research staff may have methods--mystical, time con- 


suming techniques--for coaxing approximate objective functions out 
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of reluctant decision makers. But by now the Dean is experiencing 
growing discomfort about the apparently abrupt passage from grap- 
pling with abstract and perhaps esoteric satisfaction functions 
directly to the computer and thence to "the solution." He can 
hardly help but feel that he has had scant opportunity to fully em- 
ploy his own expertise and intuition-<that, in fact, he has now 
begun to lose control over the decisions for which he is responsible. 
The Dean can avoid the trials of developing a comprehensive 
objective function end recover some of that lost initiative if he 
and his staff decide to attack this analysis in a different manner. 
Instead of formulating a mathematically complete problem and turning 
it over to the computer for solution, the Dean may deal with selected 
portions of the solution process himself, leaving the remainder to 
the machine. Computer solution techniques do not make use of all of 
the implications of an objective function anyway; loosely speaking, 
they merely interrogate it about the situations arising at particular 
alternatives. The decision maker may substitute his own judgments 
concerning particular alternatives for that part of the solution 
technique in wnich the computer would interrogate the objective 
function. The advantage of this substitution is that the objective 
function need never be specified precisely. Geoffrion, Dyer, and 
Feinberg? call this “interactive multicriterion optimization." A 
major drawback is that the type of introspection required is not as 
natural as the simple comparison of alternatives, for example, as in 


the use of an input-output model. It also involves a repetitious 





cycle of introspection and computer operation, a requirement which 
will almost surely be e or time consuming or both. Dyer* 
has implemented interactive multicriterion optimization in a form 
requiring only that the decision maker be able to choose from among 
specific alternatives. Much indoctrination of the decision maker is 
avoided thereby, at the expense of increasing the computer-decision 


maker dialogue. Saska? and Benayoun and Tergny* 


provide somewhat 
simpler interactive optimization techniques for use when the 
functions f,(x) and g;(x) 811 happen to be linear. 


Another possibility is "goal programming" ,9»9,7,8 


In using 
this approach, the decision maker sets desirable levels of accom- 
plishment ("goals") for each criterion. Linear programming tech- 
niques may be used to minimize the total weighted underachievement 
of these gcals, or to accomplish as many goals as possible in a 
specified order of importance, or some combingtion of these two. 
setting goals and ordering priorities requires some finesse; ordi- 
narily several successive modifications to the initial goal and pri- 
ority set will need to be made.9 The planning process using this 
technique is not unlike the use of an input-output model. There is, 
for example, no guarantee that the Dean will have the wisdom, 
patience, skill, or time to arrive at the most appealing possible 
projection. 

This thesis offers another way for the Dean to approach his 
planning problem. The linchpin of the method is the notion of 


“efficiency” or “Pareto optimality,” after the Italian economist who 





originated the concept.” 

Briefly, an alternative is "efficient" if it is achievable 
and if no other achievable alternative scores at least as well wih 
respect to every criterion and is actually preferred with respect to 
at least one. In tne notation presented earlier, decision set x is 


ረ 
efficient if and only if g,(x) = 0 for all constraints i, and if 


۱۱۳ 


there does not exist another decision set x' such that gi(x') 
for all i and such that f, (x') ? fi(x) for all i, and fi(x') is 
greater than f; (x) for at least one i. 

It is obvious that if the Dean is to judge his projections 
solely on the basis of their scores with respect to the criteria that 
he has laid down, then the decision set which is most satisfying must 
be an efficient decision set. The method proposed is to cull the 
efficient decision sets from the mass of all possible decision sets 
and to present them, in a clear fashion, to the decision maker for 
his perusal. ¿de is then free to consider theseralternavives in 
light of their levels of attainment with respect to the various cri- 
teria laid down. The later chapters of this thesis may offer some 
clues to now this might be accomplished, but no firm prescriptions 
ን ን ንው 20551515, for balancing imeommensureables is the focus 7 
managerial art. The signal advantage of this planning approach is 


that thas decision maker retains control over the weighing of incom- 


2 
mensur&bles in his chosen field. It is this part of the decision 


process for which hs, above all others, is suited, and for which he 


alone is responsible. 
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The disadvantages of the method presented are that the cri- 
teria and constraints must be linear, and that sometimes the set of 
efficient possibilities is not very informative. The latter case 
can arise in a number of ways. Some problems are encountered in 
which every possibility is efficient; then use of the proposed 
method does not reduce the field of possibilities from which the 
decision maker must choose. It can also happen that the set of 
efficient possibilities, though smaller than the set of all possi- 
bilities, is still too large for the decision maker to search. 
Finally, even if the set of efficient possibilities is not excessive- 
ly broad, so much detail may be required to express it that the de- 
cision maker may still not be able to cope with the flood of data he 
receives. But if one of these unfortunate eventualities should 
occur, the difficulty will be apparent before the Dean invests 
apprecisble time in weighing the possibilities. Moreover, should he 
then turn to one of the other approaches mentioned earlier, he will 
probably find that he can employ the criteria and constraint expres- 


Sions already developed. 





1-2 MODELS FOR PLANNING IN HIGHER EDUCATION - A REVIEW OF THE 
LITERATURE 

Most planning aids for institutions of higher education fall 
neatly into one of two groups. On the one hand are a host of 
special nurpose techniques addressed to problems of daily operation, 
such as classroom scheduling and library operations. The range of 
thecry employed in these tools is extensive. Schroeder! and 
McNamara? provide current surveys. On the other hand there are a 
number of broad scops institutional planning models, almost all of 
them of the straightforward input-output variety. Weathersby and 
Weinstein! have published a survey of such models, tabulated by 
their salient characteristics. These models range from simple lin- 
ear forecasting cf total enrollments and physical space requirements 
to highly involved simulations treating enrollments, faculty, finan- 
cial aspects, physical facilities, and even research, The most ad- 
vanced of these models are able to handle levels of aggregation, or 
rather disaggregation, right down to the individual lecture or 
laboratory session, and require data--masses of it--in corresponding 
detail. Apparently all models presently require their data to be 
specially prepared, but the creators of the most complex models-- 
which can not be said to be finished, but are, rather, in states of 
continuing evolution--envision that ultimately the models will 
interface with pervasive institutional management information sy3- 
tems, thus eliminating, or perhaps camouflagins, the need for 


special data coliections. 
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The two most widely known educational planning systems are 
CAMPUS, and RRFM, both mammoth simulations. CAMPUS was originated 
by Richard Judy and Jack Levinel® at the University of Toronto. Its 
latest version is the most dstailed educational planning model cur- 
rently in operation. In technical terms, it 1s simply a streight- 
forward resource-costing model, albeit a vast one. In 1970 it 
incorporated student flows, space and personnel needs, and support 
costs. Efforts have since sought to incorporate debt financing and 
to adapt the system to community colleges, state college and uni- 
versity systems, colleges with entirely individualized instruction, 
and primary and secondary education scnool systems. The model itself 
is resident on a large time-sharing network, and each terminal is 
equipped with a cathode-ray tube display in addition to normal input- 
output devices. The level of aggregation in each of the sectors 
mentioned is at the discretion of the user. In its least aggregate 
form, the model's data requirements are monumental; increasing the 
agzregation effects & corresponding reduction in the data needs. 
Fortunately, CAMPUS is designed to accommodate data in the forms in 
which it is ordinarily collected. 

RRPM, the "Resource Requirements Prediction Model",i* was 
originally developed by Mathematica for WICHE, the Western Interstate 
Commission for Higher Education. Since its inception, several ver- 
sions have been developed and the RRFM models and much related work 
have been ceded by WICHE to NCHEMS, a new "National Center for Higher 


Education Management Systems at WICHE" funded by the federal 
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government. RRM, too, is basically an ordinary resource-costing 
model. As its name implies, the principle focus of the RRA family 
is upon matters of institutional finance. RRM, however, is but one 
of a galaxy of projects planned or underway at NCHEMS; these include 
work on student and faculty dynamics, on the outputs of higher educa- 
tion, and on classification of higher education activities and 
development of precise lexicons of terms for higher education manage- 
ment. Each of these can be expscted to have impact on future ver- 
sions of RRA. The data requirements for RRPN are compatable in for- 
mat with the requirements of the Department of Health, Education, and 
Welfare’s HEGIS ("Higher Education General Information Survey”) 
report. 

Two interesting models which are not of the ordinary input- 
output variety are the goal programming model of Lee and Clayton, ® 
and the application of interactive nonlinear programming to the 
operation of an academic department, by Geoffrion, Dyer, and Fein- 
berg.* Goal programming was touched upon in the last section-- 
briefly, but in sufficient detail for present purposes. The reader 
will recall that goal prozramming is not a pure optimization tech- 
nique. It is, rather, a sequence of suboptimizations dependent upon 
the specification of the goals. Although the procedure does proceed 
to optimize once goals are set and ranked, the technique does not 
include a formal procedure for adjusting the goals settings, weight- 
ings, and rankings. Instead, these adjustments are left to the 


judgment of the decision maker, just as the decisions themselves are 
3 





وس ات س 
umn ny ee‏ 
i n itle qi 0n co ut comi‏ ےہ 2:25 መ=ጡጫመ‏ 
سیب  -‏ رم مس سسسس راس — — 
AY‏ > کم LOO 6” dm‏ 
See -“፡*፡፡መ፦››-.፣ — LO aÁ-‏ 
ae” eee OP Aap a”, aie —‏ > 1 
nnnc - m—‏ > - پس تت سینا 









12 


up to the decision maker in the case of an ordinary input-output 
model. In either case, progress toward an optimum is attendant upon 
progressive refinement of the model inputs. 

On the other hand, the method of Geoffrion, et al.,? actual- 
ly an adaptation of a well known method of solution of nonlinear 
programming problems, is a genuine optimization technique. It 
differs from the origina. all machine version in that the decision- 
maker instead of the computer is asked to interrogate the objective 
function. Thus, again, the decision maker is inserted into the 
iterative cycle. This technique will be described more fully in 


section I-3. 
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I-3 MULTICRITERION OPTIMIZATION PROBLEM ANALYSIS - A REVIEW OF TEE 
LITERATURE 


The "multicriterion optimization problem" (or "vector maxi- 


mum problem") is stated? as follows: 


Maximize F(x) 


Z 
sub esta Ora = cin 


where x is an n-vector (x3,x2,...,x,), 8nd F(x), & function of x, is 

a k-vector with scalar components (f4(x),fo(x),...,fí(x)). The 

solution to this problem is the set of all x which are efficient and 
۱ 7 7 ۱ £ E 

which satisfy the constraints g,(x) * O for i = 1,2,...,m. Solu- 

tion vector x ís efficient if it satisfies the constraints and if 

there does not exist another x' also satisfying them and having the 


properties that 


f,(x") ۱۳ os ius 1,9. ا‎ 


f(x") » fi(x) for at least one such ات‎ 
Although interest in multicriterion problems is quite recent, the 
litera bure conc=snins hese problems is already voluminous. በርኩ ዉ 
bibliography and synthesis cf work in this area, see Roy.19 Roy 
distinguished between four apprcaches to the multicritericn problem, 
the second of which is “progressive definition of preferences to- 


gether with exploration of the feasible set". In contrasting this 
4 
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approach with that which begins with the construction of a unique 
functional, he notes substantial inefficiency in the latter. The 
unique functional is constructed to be able to satisfactorily rank 
all possibilities. This is frequently unnecessary; often some 
possibilities will be obviously inappropriate--it is not necessary 
for the functional to be able to correctly rank these. For example, 
it should not be necessary to deal with inefficient solutions; a 
proper ranking of efficient possibilities will suffice. But Roy 
notes that “it is often difficult to formulize the set of efficient 
actions and work with it". The algorithms of Chapters V and VI of 
this thesis offer the means to find and specify this set when the 
f(x) and g; (x) are all linear expressions in the x,. 

Roy mentions goal programming as one method of "aggregation 
of multiple objective functions ina single function defining a com- 
plete preference order". In view of the iterative way in which the 
proper (or at any rate, final) goals and rankings are typically 
developed," goal programming could just as properly be placed in 
Roy's second category--the "progressive definition of preferences 
together with the feasible set". Lee and Clayton® have constructed 
a goal programming model for a single year university planning prob- 
lem, and Lee and Sevebick’ have applied goal programming to a dy- 
namic, multiple year municiple planning problem. Both of these 
models are entirely linear. Dyer” has linked goal prosraming with 
the more formal "interactive" exploration of the feasible rezion and 


objective function approach of Geoffrion, Dyer, and Feinberg. 
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In the latter, a proven nonlinear progremming technique is 
employed in a manner requiring only “local” information about the 
decision maker's preferences. Essentially, Geoffrion et al.l sug- 
gest the following iterative cycle. Beginning at some particular 
solution, the decision maker is asked to state tradeoffs ratios 
among the criteria evaluated at that point. This information is 
translated into a direction of improvement; the computer then deter- 
mines the criteria values achieved as the solution is varied along 
the direction chosen. These results are displayed for the decision 
maker, who selects the distance along that direction which seems to 
achieve the most desirable results. At the solution so selected, 
the decision maker states new criteria tradeoffs, and the process 
continues as before. An optimum is found when there is no direction 
of improvement. 

Roy cites two other interactive approaches, 3 4و‎ both involv- 
ing an iterative interrogation/computation cycle, and both of them 
functioning in a linear programming framework. The technique of 
Geoffrion et al.! is, of course, not restricted to the linear case. 

The characterization of efficiency in the mathematicel pro- 
gramming context is admirably summarized and extended in Philip; it 
is on this paper that much of the following analytic work is based. 
The most useful result in the paper is an extension of a theorem of 
Kuhn and Tucker providing useful necessary and sufficient conditions 
for efficiency. The theorem is that if x is on the boundary ofa 


polyhedral feasible region then x is efficient if and only if some 
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nonnegative weighting of the gradients of the constraints binding at 
x equals some weighting of the gradients of all of the objective 
criteria, where these (latter) weights are all greater than or equal 
to 1. This condition can be conveniently checked by attempting a 
linear program. Much use of this will be made in Chapters V and VI. 
Philip also presents a collection of algorithms, some due to others 
and some his own, which are addressed to the problems of deciding if 
a solution is efficient, finding an efficient solution, finding if a 
given efficient solution is unique and if not finding another, and 
finding an efficient solution which maximizes some linear objective 
function. He does not address the problem of finding and expressing 
the set of all efficient solutions. 

Finally, some interesting work is being done in the area of 
dominance structures.l? Efficiency is the weakest possible form of 
dominance; the strongest is supplied by the usual single objective 
function. Much remains to be done to enable decision maker and ana- 
lyst to deal effectively with problems for which there is readily 


available a dominancs structure somewhere between these extremes, 





17 


I-4 THE CASE FOR EFFICIENT MANIFOLD PRESENTATION 


The approach to multicriterion optimization problems which 


is suggested in this thesis consists of presenting to the decision 


maker the ' 


‘efficient manifold"--that is, the set of all efficient, 


or "Pareto optimal", possibilities. For convenience, this method 


will be styled “efficient manifold presentation", or simply "EMP", 


The case for using efficient manifold presentation rather than an 


input-output model, goal programing, or straightforward or inter- 


active nonlinear programming revolves primarily around its improved 


accommodation to the difficult circumstances of the average decision 


maker, to wit: 


a) 


b) 


Instead of offering "the answer", it provides a range of 
possibilities over which the decision maker can exercise 
his insight, taking into account unquantifiable and in- 
tangible or simply unrecognized considerations. The ab- 
sence of choice in "the answer" is coming to be recog- 
nized as a major failing in much of management science's 
mathematical modeling .“? 
It allows the decision maker to apply his managerial arts 
and his intuition in a familiar context--that is, direct- 
ly to the problem of choosing the best alternative from 
among & set of real alternatives. Most other approaches 
permit only indirect influence over the course of the 


method, for example through parameter adjustments, goals 


setting and ranking, or formulation of objective 
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functions. Any of the latter is apt to seem unnatural to 
the decision maker, and each is required before or during 
use of the computer, if one is employed, rather than 
after. 

EMP does not require the decision maker to interact with 
the computer. Besides being inherently distasteful to 
the uninitiate, interaction either requires the decision 
maker to be at a computer terminal whenever he works on 
his problem (most executives have neither the time nor 
the inclination to use computer terminals), or else it 
imposes a considerable time lag while the computations 
for each new iteration are batch processed. 

Unlike the input-output model approach, but in common 
with the other approaches, EMP fosters thought about cri- 
teria as well as constraints. This enforced discipline 
can by itself be expected to be very enlightening, even 
if EMP is not ultimately pursued. 

Most importantly, EMP, if successfully employed, can 
economize on the decision maker's most precious commodity 
--his own time. Because it requires neither iterative 
computer interaction nor difficult and time consuming 
definition of a monolithic preference function, it is 
anticipated that EMP will prove to be a most streamlined 


aid for multicriterion decisionmaking. 
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Of course, EMP is no panacea. The theory in this tnesis 
only enables one to cope with problems whose variables are continu- 
ous, whose constraints are linear, and whose criteria are linear and 
monotonically desirable. In addition, if the number of criteria is 
excessively large, or if the set of efficient possibilities is very 
complex, then the EMP method may generate more information than the 
decision maker can digest. Finally, in some problems every possi- 
bility is an efficient possibility. In such a case, the ENP method 
will result in explicit expression of the set of all possibilities 
and if there are any constraints which are nowhere binding, this 
fact will be exposed. (Such a case arises in the example problem of 


Chapter IV and VII.) So same enlightenment is possible even here. 
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1-5 A SIMPLE EXAMPLE OF EFFICIENT MANIFOLD PRESENTATION 

As a rather unrealistically simple example of the use of 
efficient manifold presentation, consider again the case of the Dean 
of a small college who is deciding upon a target steady state en- 
rollment for his school. His institution produces two types of 
graduates; liberal arts graduates (enrollment in the liberal arts 
curriculum will be represented by the variable x1); and engineers 
(enrollments in the engineering program will be x9). The annual 
number of students who graduate from each program is estimated to be 
20% of the total steady state enrollment in that curriculum. The 
Dean is not indifferent to the mix of graduates; in particular, he 
feels that a heavy preponderance of one type of graduate would be 
less desirable than &n equal number of students more evenly distrib- 
uted. Consequently he has two criteria; ihe yearly number o? liber- 
al &rts students graduated, and the yearly number of engineers 
graduated. 

Three constraints limit enrollments at the Dean's college. 
First, all students must live in the dormatories, and dormatory 
accommodations--and therefore total enrollments--are limited to 
1000. Engineering laboratory space is also limited; no more than 
550 engineers can be accommodated. Finally, the institution's annu- 
al computer budget is a firm $200,000. Usage rates at the institu- 
tion seem to average $100 per liberal arts student per year, and 


$500 per engineer per yeare 
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er Desan's problem 1s*depicted graphically in Figure I 6 
alrebraically below. The total enrollments corresponding to each 
point in the graph may be read off of the axes; the graduation rates 


are 1/5 of these. 


Criteria: Max 1/5 x, (max liberal arts graduates) 
Max 1/5 x» (max engineering graduates) 
Constraints: x +x $ 1000 (dorms) 
ራ ac 
x9 7 350 (engineering labs) 


ረ 
1002; + 500x = 200,000 (computer funds) 


A 
ود‎ 22 = O 


The three constraints are plotted in the grapn; the only enrollments 
which are achievable are those whose points lie in the hatched 
region or on its boundaries. The graph also contains two arbitrari- 
ly positioned dotted lines representing the criteria. Loosely speak- 
ing, the object is to push the horizontal one up and the vertical 
one to the right until their intersection is about to leave the 
hatched region. The point so found will be only one of a number of 
possible points with the property that further improvement of either 
criterion is impossible except at the expense of the other. 

The Dean's institutional researcher wishes to present to the 
Dean the set of all efficient enrollments. The observant reader 
will recognize thet an efficient point must lie in or on the bounda- 


ry of the hatched regicn and must have the property that there is no 





25 


other point in or on the boundary of the hatched region which scores 
at least as well as this one with respect to both criteria and 
scores better with respect to at least one of the criteria. 

These efficient points are the ones lying on the heavy lines 
in Figure II. The reader should verify that for any other point in 
the hatched region, there exists another which is at least as far 
|| 171777 add co the right, and further in one or botn dirsevions. For 
example, of the five lettered points in this grapn, only two--D and 
E--ara efficient. Point A is not efficient because other achievable 
points, including E, have greater numbers of both types of gradu- 
ates; similarly, B is "dominated" in like manner by point D. There 
are no achievable points with more of both types of graduates than 
are &chieved at C, but there are some, including D, with at least as 
many engineers &ndj more liber23l arts graduates; then C is not effi- 
cient either. 

The Dean is now able to restrict his attention to the possi- 
bilities represented by the heavy lines, rather than having to con- 
Sider all of the possibilities in the hatched region. After some 
deliberation, ne chooses the enrollments set represented by the 
point E, because he feels most comfortable with the numbers of grad- 


uates of the two types wnich it offers. He does this without con- 
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this is sufficient because to every alternative in the hatched 


region but not on the heavy lines there corresponds another on the 


rs) 


heavy lines which is clearly preferable to the first. Since point 
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CEAPTER 1 


THE SCHOOL OF HEALTH SERVICES 


11-1 RATIONALIZATION AND REALLIGNMENT OF THE NON-PHAYSICIAN HEALTH 
CARE PROFESSIONS: A NEED AND AN OPPORTUNITY 

It is common knowledge that health care delivery is a 
troubled industry today. Though it is not the purvose of this 
۴۰٢٠٠٢ to air the maladies of the medical world, certain by no 
means original observations are in order. Recent decades have seen 
the svecialization of almost all physicians; this trend has left a 
huge and growing gap in primary health care. The breach could be 
filled by non-physician intermediste-level health care profession- 
als, were they numerous enough, but unfortunately they are not. In 
fact, there 15 a serious shortage of persons with baccalaureate 
level training in health care delivery. Nurses currently comprise 
the largest proportion of such professionals and nurses are indeed 
in short supply, but the most critical shortages are in the less 
well known, more recently established and more quickly growing pro- 


fessions. It is expected that much of the future unfilled demand 


“The breadth of this family of health care professionals is 
11010016 00006 0۰0710010116 list of curricula for 208ھ‎ School of 
Health Services of the University of Nevada. 
Anatomy Medicine 
Biochemistry Nursing 
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for baccalaureate level health care professionals will be for mem- 
bers of these new fields, and for professionals whose professions 
have not yet been defined. Unfortunately, the newest professions 
are unsupported by educational establishments equal to the task of 
training sufficient numbers of graduates in the new fields, and even 
the long established schools of nursing are beset both by financial 
problems and by a crumbling of the general consensus on what nursing 
education should encompass. Educational institutions which might be 
interested in contributing to the solution are handicapped by a lack 
of agreement on the character of the new professions and on the 
forms of the curricula that students of these professions should 
follow. Compounding this, the new depression in higher education, 
with its general decreases in state and federal funding support, 
discourages experimental ventures. 

ል complicating factor in the shortage of intermediate level 
health care professionals is the absence of well defined paths for 
advancement. Provision for career ladders for these professionals 
is a very difficult problem, one which will require consideraole 


study before appropriate solutions are found. Unfortunately, 


Biomedical Engineering Pathology 
Biomedical Instrumentation Physical Therapy 
Dietetics Physiolocy 
Health Education Science Writing 
Medical Technology Speech Pathology 


Even this extensive list is not exhaustive. 
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solutions to this problem and to the problems of curricula organiza- 
tion are intimately related. 

New curricula, facilities, and agreement on professional 
boundaries are needed to alleviate this manpower shortage. Because 
of the decentralized nature of higher education, these ars apt to 
develop slowly. Progress will be mcst rapid if the more prestigeous 
institutions of medical education will take the lead. As a recog- 
nized pace setter in medical education, it is fitting that the Johns 
Hopkins University should pioneer in these advances. It is the 
recognition of these needs and this mandate which led this Universi- 


ty to create the School of Health Services. 





لت 


11-2 GENESIS OF THE SCHOOL OF HEALTH SERVICE 

The Jonns Hopkins University and Hospital have long partici- 
pated in education for the allied health professions, both through a 
distinguished diploma program in nursing and through numerous ad hoc 
programs for training intermediate and lower level health care per- 
sonnel. Naturally, as in all institutions of higher education, the 
orzanization and content of curricula have been under constant study 
and discussion by a host of committees and individuals. The chain 
of events leading to the establishment of the School of Health 
Services, heiir to this long succession of interest and participation 
in non-physician medical education, may be said to have formally 
begun in 1967 when an ad hoc committee of the Faculty of Medicine 
under the chairmanship of Dr. Russell Morgan recommended the forma- 
tion of a College of Allisd Medical Sciences. This recommendation 
was approved by the Advisory Board of the School of Medicines and in 
1968 the means cf implementation became the subject of study of a 
University-wide ccrmittee chaired by Dr. Morgan. As a member of 
this committee, Dr. John P., Youngs, Associate Provost for Planning 
and Advanced Policy Studies, prepared an extensive and detailed 
feasibility study. The study, which was completed in December of 
1968, proposed the initial establishment of two baccalaureate pro- 
grams, one in nursing sciences and another in clinical laboratory 
sciences, and included estimated costs and curricula outlines. This 
report was endorsed by the Morzan committee and provided tne focus 
for further discussion and deliberation, in 1969, in two informal 


joint meetings of the Advisory Boards of Medicine and Hygiene, and 
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the Academic Council of tne Faculty of Arts and Sciences, during 
which the nascent division was given approval in principle. In 
Decerber of 1970, the Center for Allied Health Careers was estab- 
lished under the direction of Dr. Dennis G. Carlson and charged with 
continuing to plan for the education of lower and middle level 
health care professionals. The Center produced A Planning Report 
for Education and Training in Health Services in May of 1971 (re- 
vised in July of the same year) which again called for the establish- 
ment of a distinct scnool within the University, under the title 
“School of Health Services". In another significant development 
that summer the Johns Hopkins Hospital, whicn is organizationally 
distinct from the University, commenced the closure of its three 
year diploma program in nursing by not admitting students for the 
class of 1974. In October of 1971 the Board of Trustees of the Uni- 
versity approved the formation of the School of Health Services and 
on 1 March 1972, Dr. Malcolm L. Peterson was appointed Dean of the 
new school. In February of 1973 the Board of Trustees authorized 
the appointment of a faculty and admission of students for the 


following academic year. 
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11-3 CHARACTERISTICS AND CURRICULA OF THE SCHOOL OF HEALTH SERVICES 

The new School o2 Health Services fell heir to several 
health care personnel education programs which existed prior to the 
School's formation. These included an associate degree program for 
Health Assistants, operated in collaboration with Essex Community 
College, and a course of instruction for Pediatric Nurse Practition- 
ers, operated in collaboration with the Department of Pediatrics of 
the School of Medicine. 

The first new program to be initiated by the new School will 
be the Health Associates program, which will prepare individuals to 
deliver health care in ambulatory settings. The curriculum is to be 
two years in length; incoming students will be required to have com- 
pleted at least two years of college before matriculating. At the 
end of the course, a baccalaureate degree will be awarded. The cur- 
riculum for this program will be presented in a uniqus format. 
Instead of courses organized by traditional disciplines (biology, 
chemistry, etc.) material will be presented in "modules" directed 
toward commonly encountered health care problems. The typical 
module will contain elements of instruction in several branches of 
theory and in pertinent tasks and skills. Faculty who will be in- 
volved in the Health Associates program can be divided into three 
broad categories. A "primary teaching team" will be resvonsible for 
most of the instruction of each group or "section".of 25. students, 
The members of the primary teaching team will be one physician, two 


"social generalists" (anthropologists, mental health counselors, 
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social workers, health educator/community organizers, etc.), and 
three "health practitioners" (persons with approximately the quali- 
fications of the graduates of this or the nurse practitioner program 
to be described later). Supplementing the primary teaching team 
will be faculty specialists (social scientists, surgeons, patholo- 
gists, etc.) and clinical supervisors (who will oversee practical 
training on location in various health care delivery settings). 
Finally, the faculty will be supported by an extensive media re- 
sources program whose function is to amplify their effectiveness 
through the employment of modern technical educational aids. 

The second baccalaureate program planned is the Nurse Prac- 
titioner program. It is expected that this program will be very 
Similar to the Health Associates program except that its focus will 
be on health care delivery in institutional care settings rather 
than in ambulatory settings. Other baccalaureate programs under 
consideration for the immediate future are programs in Health Ser- 
vices Managenent, Environmental Hygiens, and Clinical Laboratory 
Science. 

Additional new programs are planned for the more distant 
future. It is anticipated that an undergraduate program in mid- 
wifery and an Adult Nurse Practitioner program similar to the ongo- 
ing Pediatric Nurse Practitioner program may be initiated. Also 
planned are a program in Nuclear Medicine Technology, and one in 
Clinical Laboratory Technology to be conducted jointly with Essex 


Community College in the same manner that the Health Assistant 
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program is now 6 ۰ 

While much study of the role, mission, and feasibility of 
the new School was done prior to its establishment, virtually all of 
the decisions relating to program and curricula, staffing and admis- 
sions, funding, and accession and allocation of physical facilities 
were left to the staff of the School. The urgent need for the 
School's graduates and the undesirability of a long preparatory 
period before tuition and operating subsidy grants could be realized 
dictated that the School commence operations as soon as practicable; 
accordingly, it was decided at the onset to focus on development of 
the curriculum for the Health Associates program, the first new pro- 
gram to be instituted. Because both the Health Associate concept 
and the modular teaching format are innovations, development of cur- 
riculum is necessarily slow and laborious. In addition to this cur- 
riculum development, the administrators of the School must locate 
faculty and funding, decide upon and implement the mechanics of ad- 
missions, and acquire facilities and equipment. All of these tasks 
will continue as long as new programs are being introduced (at the 
rate of one baccalaureate program annually, at least through 1977). 
Concurrent with all of these immediate responsibilities is the less 
urgent but nevertheless important problem of longer range planninz-- 
particularly the planning of relative emphasis on curricula and of 
projected requirements for personnel, facilities, and funding. 
Under the cimeumstances, this plamine- is exceptionally difficult to 


carry out. Unlike an ongoing institution, the School of Health 
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Services has no historical basis from which to extrapolate into the 
future, no current faculty or students to influence ranges of 
choice, no traditional commitments for facility or financial sup- 
port. Moreover, because the very curricula are new and the mode of 
instruction is unprecedented, only the most general sorts of infer- 
ences can be drawn from the experiences of other institutions. Also, 
curricula formalize and specify the relationships between students 
and staff, so that at least this important information will not be 
known with certainty until all curricula are developed--an eventu- 
ality which is years in the future. Nevertheless, decisions will be 
made in the interim, if only by default, and they can only be made 
with the imperfect information available. 

One set of longer range decisions to be made involves 
enrollments and staffing levels. In particular, the distribution of 
the School's limited resources over the various baccalaureats pro- 
grams must be made. A number of planning approaches potentially 
useful in higher education were listed in Chapter I; of these, 
efficient manifold presentation is perhaps the one best suited to 
these particular circumstances, for the following reasons. To begin 
with, the relatively small size of the School makes complex computer- 
ized input-output models a poor choice. Since the staff of the 
School is working under the press of many immediate problems, the 
computer-interactive approaches may be ruled out as too time consum- 
ing. Similarly, there is scant opportunity for the difficult 


development of a composite objective function for use in ordinary 
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nonlinear programming optimization. Efficient manifold presentation, 
on the other hand, suffers none of these drawbacks. It has tne 
potential to provide significantly better insight than scratch md 
input-output projections for very little additional decision maker 
effort. 

Tne next chapter contains an assortment of components suit- 
able for constructing linear multicriterion educational planning 
models--that is, planning models for higher education which are 
amenable to efficient manifold presentation. The following chapter 
contains a model tailored to the circumstances of the School of 
Health Services. The solution technique is developed in Chapters V 
and VI, and the results of its application to the School of Health 


Services model are given in Chapter VII. 


- 
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FOOTNOTES 


Ismith, G. T. and N. A. Baxter, University of Nevada Medical 
Education Feasibility Study, Nevada State Bureau of Business and 
Economic Research, November 1968. 
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CHAPTER III 


A FRAMEWORK FOR ENROLLMENT AND FACULTY PLANNING 


111-1 GENERAL 

This chapter contains a general model of certain aspects of 
the operation of an educational institution and a broad brush 
characterization of the use of the "efficient manifold presentation 
method" to augment the intuition when considering admissions and 
staffing decisions. In Chapter IV, this model will be tailored to 
mit the circumstances of the School of Health Services. As an 
illustretion of the use of the efficient manifold presentation 
method, the results of its application to the model of Chapter IV 
are presented in Chapter VII. 

It is natural to think of planning problems in terms of (1) 
ends or goals, and (2) circumstances prescribing the allowable 
actions which could be employed in the pursuit of these goals. The 
reader will recall from Chapter I that the actions of interest here 
are student admissions and faculty staffing decisions. These, the 
"variables" for the following model, will bəs represented by the syn- 
bols c and f, appropriately subscripted. The goals and circum- 
stances affecting or affected by these variables are the criteria 
and constraints. In the following two sections, some criteria and 
constraints wnich are likely to arise in admissions and staffine 


planning for institutions of higher education are discussed. 
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III-2 SAMPLE CRITERIA 
From the, modeler's point of view it is quits unfortunate 
that of the many advertised goals of the world of higher education, 
few seem to be susceptible to precise definition or accurate measure- 
ment. Nevertheless, the following list of candidate criteria is 
offered. This does not pretend to be a comprehensive list of all 
the motives that propel educators; it doss purport to be a reason- 
ably complete list of the measurable goals whose achievements are 
noticeably affected by decisions on numbers of admissions and 
exmeulty. The criteria are as follows: 
Graduate as many high quality students from a given program 
as possible (applicable to each degree program) 
Maximize the excess of income over expenditures--or minimize 
the shortfall 
Minimize the size of the average class in a given category 
of courses (applicable to each course category) 
Minimize the deviation fran ۵ specified ideal discipline 
distribution within the faculty 
Minimize the deviation from a specified ideal curricular 
distribution within the student body 
Minimize faculty terminations 
Minimize the additional capital investment required to sup- 
port the chosen enrollment-staffins plan 
Pemoresdtacussing each of these 117 


Sane general remarks are in order. First, none of these criteria is 
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intended to represent, by itself, a capsule pnilosophy of management 
for any educator. It must be borne in mind that these criteria are 
members of a family of yardsticks with which tne thoughtful adminis- 
trator measures the performanca or anticipated performance of his, 
and other, institutions. None of these criteria claimsto provide 
any ordering relationships by we except wnen all other cri- 
teria are equal. As the explanations of these criteria and of the 
constraints in the next section unfold, the opportunity will be 
taken to introduce equations and some occasional notation, in antici- 
pation of the formal mathematical model of Chapter IV. Some of the 
criteria above are stated as quantities to be maximized; others as 
quantities to be mininized. For uniformity, and to be consistent 
with the theory developed later, all of tne mathematical expressions 


for these criteria will be expressed as quantities to be maximized. 


Graduation rates: 

Heading the list of criteria is that of graduating as many 
high quality students as possible from a given program. The virtue 
of this--wnich used to be self-evident--is somewhat tarnished by tne 
present surfeit of graduates in a number of fields and by the peren- 
Pla I2 few; nevertheless this 1s a plausiole criserion if 
only because the just, equitable, and economical] way to limit output 
is to put additional emphasis on other worthy and partially conflict- 
ing criteria, sucn as annual surplus. Additionally, when all other 


things are equal, it is almost certainly the case that in the eyes 
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of any particular institution, more of its own graduates is pre- 
ferred to fewer, even when the market is flooded. 

The mathematical quantity to maximize is the product of 
enrollment in the program, or if enrollments are categorized by pro- 
gram year then the enrollment in the last year of the program, times 
a constant representing the proportion of the enrollment who become 


graduates each year. 


Annual financial surplus: 

In the literature of institutional research, one occasionally 
finds the university described not as a "non-profit" institution but 
as a "not-for-profit" institution. The distinction is more important 
and less subtle than it appears. The days of profit oriented pro- 
prietary schools are indeed gone, but the need for fiscal responsi- 
bility is not. Financial surplusses are still desirable, for they 
can be devoted to any number of desirable uses, from improving 
quality to increasing endowment. The second criterion is to maxi- 
mize profits; this enables one to distinguish between enrollment- 
staffing plans which are indistinguishable except with regard to 
expense. The mathematical quantity to be maximized is the sum over 
all students of the net income per student less the sum over all 
faculty of the net expense per faculty member less a constant sub- 
suming all expenses and incomes not proportional to students or 
faculty. 


Maximizing the sum excluding this last constant amounts to 
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the seme thing. Significance is to be found in the relative rather 


than absolute valuss of the scores according to this criterion. 


Student /faculty ratios: 





The third criterion is to minimize average class sizes. The 
belief that knowledge is better transmitted to small groups than to 
large groups is à prominent feature of educational lore; it is 
expected that when all other things are equal, enrollment-staffing 
plans with small classes are preferred. The natural way to express 
the relationsnip between students, faculty, and class sizes would be 


the following: 


average class size = (constant)x(numbers of students) 


(numbers of faculty) 


where the constant subsumes average course loads and faculty teach- 
ing loads and the like. Unfortunately this equation is not linear 
in the two variables, and linearity is one property which is essen- 
tial for the analysis of Chapter ۷٢. One way to deal with this prob- 
lem is to use two separate criteria: numbers of faculty (to be 
maximized), and the negative of the numbers of students (to be maxi- 


mized). 


Faculty composition: 
The composition of the faculty by disciplines is important 


for a number of reasons. The most obvious is that the faculty in 
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each discipline should be numerous enough to teach the courses in 
that discipline which are needed to meet the requirements of the 
students enrolled. This requirement will be expressed as a con- 
straint in the next section, Another family of constraints could 
arise fran considerations, such as accreditation standards, requir- 
ing at least a certain proportion of the faculty to hold particular 
types of degrees. In addition to these, the disciplinary composi- 
tion of the faculty is important in its own right. For example, 
virtually every institution will insist on maintaining a humanities 
program even if it is not required by one of the constraints just 
mentioned. The reason is the humanizing effect these faculty are 
expscted to have on students and on their own more technical col- 
leagues. If the Dean is able to express his notion of the ideal 
community of scholars in terms of percentage composition, then the 
deviation fran that ideal can be expressed and minimized. Suppose 
that it is desired that each faculty type should comprise the pro- 
portion b; of the entire faculty. Introduce the following con- 


seraints, one for each category 1: 


- c 
b. 7, f;-fitxpoi > 0 
i 


where f; is the number of faculty of category i, and xi and x; are 
the positive and negative differences between the actual and desired 


proportions, and where all of these variables sare constrained to be 


nonnegative. Then the quantity to be maximized is; 





Maximize 22 -x.-Xj 
i 
Involuntary faculty terminations: 

One of the most difficult, painful, and controversial 
activities which circumstances sometimes force a Dean to carry out 
is the task of reducing the size of the faculty. A most reasonable 
objective is to reduce as much as possible the expected number of 
involuntary (nonretirement) terminations which are built into any 
given enrollment and staffing plan. The number of faculty of cate- 
gory i in a given year k can be approximated by the following equa- 


tion: 
5 
44 ጄ “ዶሪ , دم وگ‎ hike Mik 
T 


where fi, is the number of faculty of type i in year x, A is the 
proportion of a previous year's faculty wno remain voluntarily or 
wish to do so for the year in question, hi, is the number hired, and 
riy is the number whose terminations begin with the year k. If this 
constraint is added for each faculty category, then the quantity to 


be maximized is: 


Maximize 2. “Iai 
1 


Capital investment: 


The last criterion on the list is to minimize the amount of 
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additional capital investment (e.g. new construction) required. If 
Of and Ps, are the amounts of facility type k required per faculty 
member of category 1 and per student of variety j respectively; if 

Ex is the amount of this facility already available; and if dy, is 


the capital investment per additional unit of this facility, then 


adding the following constraint for each k 


2.06 ٢ 2 n C 47Xk = By 


(x, ue) 
allows uss of the following criterion: 


Maximize = “Ay Xx 
k 
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111-3 SAMPLE CONSTRAINTS 


A few constraints have already been mentioned in conjunction 


with the criteria of the last section. Some additional constraints 


ares 


Enrollments, numbers of faculty, and student-faculty ratios 
are related through a mathematical identity. 

Admissions and enrollments are matnematically related. 

The number of applications sets a limit to the number who 
ultimately matriculate» 

The number of prospective faculty available limits the 
amount of hiring which may be done. 

The facilities availabe limat enrollments amc ከ. 

Upper (lower) bounds on the proportion of the faculty satis- 
fying certain characteristics may be imposed. 

Upper (lower) limits on the proportion of enrollments con- 


posed of certain varieties of student may be imposed. 


Student/faculty proportionality: 


The first two constraint types are actually mathematical 


identities expressing obvious but occasionally necessary relation- 


ships. 


In the last section it was noted that student-faculty ratios 


are of great importance, and a pair of criteria based on minimizing 


the class size was proposed. In the event that it is desired to fix 


the class size, tnen the following constraint is in order: 
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Jfaculty = (#students)x(rate of course registrations per student) 
class stu-fac ratio)x(üfcourses per professor 


The only variables here are numbers of students and faculty, so this 
is linear. When faculty is classified according to discipline and 
students according to curriculum or according to curriculum and year 
of curriculum then the last equation becomes a family of equations, 
one for each faculty category. The matrix of course registrations 
by discipline and curriculum, or by curriculum and year of curricu- 
lum, is called the "induced course load matrix". This "ICLM" is an 
almos?; ubiquitous feature of higher education planning models; it is 
frequently the Achilles heel of the model, too, because when the 
curricula have room for electives, the entries in the matrix devend 
on the mercurial interests of the student body. That is, when the 
Ser ricula are flexible, the ICL is difficult to establish and prone 


to surreptitious change. 


Relationsnips among matriculations and enrollments: 
AN tn ی اس ار‎ ዘ شا نتا‎ 


The second constraint may be expressed: 


mumber of 
students - (number this year) x (proportion + (new students) 
next year who stay 


This constraint becomes a family of constraints when the student 
body is categorized. The proportionality figures here will be know- 


able to a greater or lesser extent depending on whether curricula 
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guidelines are firm or flexible and whether it is difficult or easy 


for the student to transfer from one curriculum to another. 


Staff and student availability: 

The next two constraints say simply that no more students 
(faculty) may be matriculated (hired) than are available. Some 
authorities forecast a decline in college attendance in the coming 
decade, so matriculations constraints are apt to be realities for 
many institutions. Despite the very visible abundance of advanced 
degrees in some fields, genuine faculty availability limits will not 
be unknown either. Such a case arises in the example of the next 
chapter. If psy is the number of type i students matriculated in 
year k and Piy is the number of bona fide applicants who ars accept- 
able that year, then the expression is 


4 
Pik 7 fik 


The analogous faculty inequality is 


c 


hix 7 Huy 


where hs; is the number of category i faculty hired in year k and 


H., is the number of such faculty available then. 
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Facilities: 

Associated with any operation of the size and complexity of 
an institution of higher education will be myriads of physical 
restrictions that the operation must stay within. Most of these 
never prevs a hindrance because other restrictions take effect first 
or because the desired course doesn't happen to pass anywhere near 
this particular shoal--but it is not always possible to predict 
which restrictions are going to prove to be real, and which are not, 
so it is advisable to add the not unlikely ones to the model just in 
cases Scme of the more likely to be encountered are numbers of 
classrooms, of laboratories, and of offices. Medical curricula have 
special requirements and, correspondingly, special constraints such 
as limited numbers of patients, wards, and examining rooms. The 


standard expression for such a constraint is: 


د 
avaiiability = sum over all users of each user's demand for the‏ 


۱ 


Upper and lower bounds on types of faculty: 
Bounds on certain components of the faculty might arise, for 
exemple, from accreditation requirements. Whatever the motivation, 


the appropriate expression for an upper bound is: 
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where g; is the proportion of faculty of category i who have the 
characteristic of interest, and G is the specified upper bound on 
the proportion of the entire faculty who share this characteristic. 
In the case of a lower bound, the sense of the inequality is 


reversed. 


Upper and Lower bounds on the curricular composition of the student 
body: 
This last type of constraint is of exactly the same form as 


that just discussed. An upper bound constraint is: 
ae 
ረ, (94-9)0; 0 


where 54 is the enrollment in curriculum i, and q and Q refer to 
students in a manner analagous to g and G above. For a lower bound, 


the following constraint is in order: 


S 
». (94-9)0ፈ ። 0 
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111-4 A CHARACTERIZATION OF THE EFFICIENT MANIFOLD 

The uses and limitations of the "efficient manifold presenta- 
tion" method were summarized in Chapter I, section 4. In this 
section, the form of the efficient manifold will be delineated, a 
method for finding and characterizing it described, and the rela- 
tion between the form of the efficient manifold and the character of 
the model from which it derives will be discussed. Some of the lan- 
guage of linear programming will be employed, but rigor is postponed 
until Chapter V. 

Imagine that the Dean and his staff have formed a model of 
their institution, cobbled out of the raw materials presented in 
sections 2 and 3 of this chapter, with perhaps some additional cri- 
teria or constraints to fit the particular circumstances of their 
institution. This model is entirely linear. The question they must 
now address is how to find the set of all efficient solutions, the 
"efficient manifold." 

Recall the definition of an efficient solution; it is a solu- 
tion to the constraint equations which has the property that there 
exists no other solution to those equations (in the linear program- 
ming vernacular, no other "feasible" solution) which is at least as 
desirable as this one with respect to each separate criterion and 
more desirable according to at least one. Now suppose that all of 
the criteria in the model are added together to form one objective, 
and let us suppose that the linear program with this objective and 


with the constraints of the Dean's model is solved. If there are 
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solutions but no finite optima, then the model is badly formulated, 
for no real world problem with physical interpretations can have an 
unlimited optima. If there are no solutions to the constraint set, 
then either the model is poorly formulated or the situation of the 
institution is hopeless. If there is a finite optimum, then this 
optimal solution is efficient. For if it were not, there would 
exist another feasible solution with the property that it scores at 
least as well with respect to each criterion and better with respect 
to at least one. Adding these scores together gives a score for the 
linear program objective which is greater than that of the supposed 
optima, which is of course impossible. 

This first efficient solution happens to be a vertex of the 
feasible region--a "corner" on the polyhedron of possiole &lterna- 
tives, because it is a basic solution to the set of constraint equa- 
tions e 

A consequence of & theorem in Chapter V is that either the 
efficient manifold comprises the entire feasible region or else the 
efficient manifold is contained entirely in the bounding hyperplanes 
(surfaces of any dimension less than the dimension of the problem 
space) of the feasible region. More than that, the efficient mani- 
fold can be expressed as a collection of polyhedra formed by the 
intersection of the feasible region with supporting hyperplanes 
(thus if the feasible region is three dimensional, its efficient 
manifold--presuning that it is not the entire feasible region--is 


some collection of the faces, edges, and vertices of the feasible 
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region). Incidently, every corner of such a polyhedron formed by 
the intersection of ths feasible region and a supporting hyperplane 
is also a corner of the feasible region; thus the intersection can 
be represented as the set of all convex combinations of some set of 
efficient vertices of the feasible region. 

Another consequence of the theory of Chapter V is that the 
efficient manifold for such a problem must be connected. From this 
it is clear that any efficient vertex may be reached from any other 
efficient vertex, such as the one found by ths linear program above, 
by passing along seme chain of neighboring efficient vertices 
(neighboring vertices are vertices connected by an edge). Then the 
way to find all efficient vertices is clear; simply check for effi- 
cency the neighbors of each efficient vertex found. It happens that 
it is a simple matter to determine whether or not a vertex is effi- 
cient. 

Having found all of the efficient vertices of the feasible 
region, the next step is to use these to express the efficient mani- 
fold. It has been noted that each efficient polyhedron in the effi- 
cient manifold can be represented by a set of efficient vertices. 
Then the efficient manifold can be represented by one or more sets 
of efficient vertices, where each set represents all convex combina- 
tions of the members of the set. Not all sets of efficient vertices 
specify efficient polyhedra, but a theorem attributed by Philip! to 
Kuhn and Tucker provides the ability to recognize whether or not & 


given set does. This theorem also provides the means for 
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ascertaining the efficiency of a given vertex. One procedure for 
searching through the possible combinations of efficient vertices is 
presented in Chapter V, and several more are presented in Chapter VI. 

One rather troublesame difficulty which may arise is that 
the efficient manifold and the feasible region may turn out to bs 
identical. It will be show that a sufficient condition for this to 
be the case is for the set of criteria to be such that some positive 
weighting of them is identically zero. In this case, any arbitrary 
solution in the feasible region must be efficient because motion 
away in any direction is "down hill” with respect to at least one 
criterion, or if & direction is found which is not down hill with 
respect to any, then it is not up hill with respect to any either. 
This condition is easily checked. 

Naturally, the efficient manifold is sensitive to the skill 
with which the model is: built. So long as the criteria and con- 
straints are correct and complete, the true optima will be in the 
efficient manifold somewhere. As in all mathematical programming, a 
poorly constructed constraint set may exclude the true optima, or it 
may include too much and cause the decision maker to settle upon an 
unattainable solution. If a criterion is omitted, the efficient 
manifold may be smaller than would otherwise be the case, and the 
true optima could be omitted. If too many criteria are employad, 
the efficient manifold may turn out to be unnecessarily large, and 
in unusual circumstances it may also eliminate portions of the effi- 


cient manifold arising from the problem with fewer criteria. In 
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this latter event, the new efficient manifold will still contain 


efficient solutions with scores identical to the omitted ones, 
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FOOTNOTES 


Philip, J. “Algorithms for the Vector Maximization Problem", 
Mathematical Programming, III (1972), pp. 207-229. 
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CHAPTER IV 


THE SCHOOL OF HEALTH SERVICES EQUILIBRIUM PLANNING MODEL 


IV-1 GENERAL 

The previous chapter outlined an assortment of linear cri- 
teria and constraints which might prove useful in modeling enroll- 
ment and staffing decisions for an educational institution. In this 
chapter, some of these are used to build a model of one particular 
institution, the newly established School of Health Services of the 
Johns Hopkins University. Some of the history and characteristics 
of this interesting educational venture were summarized in Chapter 
II. The model herein was prepared primarily in order to provide a 
test of the analysis technique suggested in this thesis., It must be 
understood that although the author received considərable assistance 
from the Dean and members of his staff, neither the model nor its 
data bear any University endorsement, official or otherwise. 

Several conceptual models of the School of Health Services 
were formulated before the following equilibrium model was decided 
upon. This model has several prominent advantages. It is mathe- 
matically much gnaaller and less complex than a dynamic model, and 
being a steady state model, it can not suffer transient effects. 


Perhaps most importantly, its outputs are much briefer and easier to 
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comprehend than are those wnich a dynamic model would yield. On the 
negative side, steady state is a prospect but dimly perceived by the 
staff of the Scnool; growth is the focus of attention. Nevertheless, 
a steady state projection has some significance. It can expose the 
limits to growth which are posed by the resources available to the 
School and reveal which resource limitations will indeed operate to 
deter further growth. Knowing this, the Dean would then be enabled 
to direct his attention to ameliorating the potentially confining 
circumstances in advance. Financial projections for these steady 
states will also be illuminating. Finally, the disciplined and sys- 
tematic enquiry which goes into developing such a model is valuable 
on its own account. 

Because this model is a steady state model, some simplifica- 
tions are possible over the dynamic case, which is implicit in some 
of tne expressions and equations of the previous chapter. A steady 
state nature is imparted by presuming that faculty do not age or 
retire, and that student matriculations are just sufficient to bal- 
ance graduations and dropouts, so that the total enrollments, and 
the enrollments in each year of each curriculum, are the same from 
year to year. The latter presumption permits the representative 
student's entire career to be telescoped into a single entity. Thus 
the model to follow counts only beginning students in each curricu- 
lum; each such student is credited with the tuition income he will 
generate and the resources he will require over the course of his 


passage through the School. The model is further simplified in that 
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student related variables and constants need only be differentiated 
according to curriculum, rather than by curriculum, by year within 
the curriculum, and by calendar year. 

The model will be presented in three segments. General 
remarks and assumptions will come first; then criteria will be pre- 
sented. Thirdly, the constraints are discussed. Data will be pre- 
sented where it first arises incident to the criteria and con- 
straints. After the discussion of the constraints is concluded, the 


mathematical formulae of the model will be summarized. 
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IV-2 ASSUMPTIONS AND NOTATION 

On the advice of the staff of the School, the scope of the 
School of Health Services equilibrium planning model is limited to 
consideration of enrollments and staffing for the Health Assistant 
program and the impending baccalaureate programs for Health Associ- 
ates, Nurse Practitioners, students of Health Services Management, 
and students of Environmental Hygiene. Each of these curricula is 
or will be two years long. 

One very expensive class of faculty does not appsar anywhere 
in the model to follow. These are the "faculty specialists” --the 
surgeons, microbiologists, pharmacologists, etc. They are excluded 
because in the judgment of the staff of the Scnool their numbers 
will not vary significantly with changes in enrollments. Because of 
their high price and because instruction is not their only function 
(another primary responsibility is curricula guidance) they will bs 
staffed in sufficient variety to represent the relevant specialties 
and in sufficient strength to present, without duplication, the 
relevant material in their specialties. Every effort will be made 
to adjust to changing enrollment patterns by changing the sizes of 
the groups which they address. 

It was noted in Chapter II that primary teaching teams will 
be responsible for groups of 25 students each. There is one excep- 
tion--the Health Assistant program. In this program, each primary 
teaching team will be responsible for a thirty student "section", 


On consultation with the staff of the Scnool, it was decided 
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to presume that no students will drop out, either temporarily or 
permanently. It is believed that careful selection of applicants 
and close counseling of enrollees will accomplish this. It is also 
presumed, on the strength of the diversity and relatively short 
duration of these programs, that--with one exception--there will be 
no student transfers among programs. The exception arises in the 
case of the Health Associate and Nurse Practitioner prozrams. These 
two programs differ principally in their emphases--the one being 
primarily concerned with ambulatory care, the other with institu- 
tional care. Current data estimates for the two programs are iden- 
tical because the staff believes that the differences in emphasis 
will not be reflected in costs and facilities usages. Consequently 
these programs would be indistinguishable in the model, so their 
matriculations variables have been merged. According to the 
School's staff, a merger in fact is a distinct possibility. 

The variables used in the School of Health Services Equi- 
librium Planning Model are as follows: 

C, = Annual matriculations in the Health Assistant program 

Co = Annual matriculations in the Health Associate and Nurse 

Practitioner programs 
Cz = Annual matriculations in the Health Services Management 
progran 
C4 = Annual matriculations in the Environmental hygiene pro- 
gram 


Y4 7" The number of offices converted to dry laboratories 
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Y, = The number of hours weekly during which large class- 
rooms are utilized in a partitioned mode (that is, used 
with the room divider closed, for two different 
classes). 

Yz = The number of hours weekly during which large class- 
rooms are filled by (single) medium or small classes.» 

Y, ^ The number of hours weekly during which medium class- 
rooms are filled by small classes. 

The symbols a and A, appropriately subscripted, will repre- 
sent constants on the left and right sides of the constraints. The 
first subscript beneath each is the number of the constraint; that 
is, all a's and A's with the same first subscript belong to the same 
constraint. Where a second subscript occurs, its significance will 
be apparent from the accompanying text. The constants T; represent 
income less expenses associated with matriculated students; the 
derivation of the T; from familiar data is expleined in section 
IV-5, in conjunction with the annual budget surplus criterion. 

No faculty variables are required because no faculty cate- 
gories are anywhere substitutable for one another. Consequently, 
given enrollments may be translated into exact faculty requirements, 
and vice versa. The known relations between enrollments and faculty 
are used to attribute faculty related requirements directly to the 


student variables. 
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IV-3 CRITERIA 

Two of the sample criteria listed in section III-2 were 
selected for use in this model; they are numbers of graduates of 
each curriculum, and annual budget surplus. The remaining candi- 
dates are inappropriate for a variety of reasons. Faculty termina- 
tions can be projected only from dynamic models; this one is static. 
The unique modular method of instruction contains within it 
decisions on student/faculty ratios and on faculty discipline dis- 
tributions, so these, also, are inappropriate criteria. Lastly, the 
55511 of the School feels that there is little near term prospect of 
a building program expressly for the School; the capital investment 


criterion is thus superfluous. 


Numbers of graduates: 

Ali otner things being equal, more graduates of a given cur- 
riculum are preferred to fewer. Since it is assumed that there will 
be no dropouts, the steady state number of graduates per year is 
equal to the number of students matriculated per year. Thus the cri- 


teris are: 


Maximize Bealth Assistant program gracuates = Max Ci 
Maximize graduates of the Health Associate and 
Nurse Practitioner programs = Max Co 


Maximize graduates of the Health Services 


ዘ 


Management program Max Cs 
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Maximize graduates of the Environmental Hygiene 


program = Max Ca 


Annual budget surplus: 

Because net financial position is one of the major concerns 
associated with the birth of a new enterprise, the criterion "maxi- 
mize instructional revenues less instructional expenses" is chosen 
as one of the criteria for this model. The most desirable financial 
measure would be the net annual financial position of the School, 
but, as noted in Chapter III, the sum of tuition and other student 


related incomes less the ' 


"proportional" costs (those which vary in 
proportion to numbers of enrollees) of instruction will suffice. 
The natural expression for this quantity would be the sum of the 
individual contributions of each student enrolled and each faculty 
member present during a given steady state year. But enrollments 
and numbers of faculty are not variables of this model; only 
matriculations are readily available. However, total enrollments 
and faculty size are directly related to matriculations, so the 
figure desired can be had by attributing the costs to the enroll- 
ments and faculty as before, and then attributing these to matricu- 
lations. Since it is assumed that there will be no dropouts, cor- 
responding to each newly matriculated student will be one enrolled 
second year student. An equivalent approach would be to calculate 


the costs and benefits for each student throughout his passage 


through his curriculum; in this approach, the first and second years 
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are thus represented just as they would be if the one year contribu- 
tions of the first and second year students simultaneously present 
in a given year were computed. This device of attributing to the 
matriculating student his entire curricular requirements will be 
employed several times in this chapter. In this instance, it will 
of course include those figures pertaining to faculty whose presence 
is contingent on the student's presence somewhere in his curriculum. 


The criterion is: 


4 
Maximize 2, TC, 
i=l 


where T; is the income less expenses associated with each student 


matriculated in curriculum i. These are estimated to be: 


T4 = 80 

To = $-11,962 
Ts = $-4644 
Tg = $-4734 


The first is presumed to be zero because it is expected that the 

National Institutes of Health will continue to fully fund the Health 
Assistant program. The other coefficients contain tuition and vari- 
able costs which are directly attributable to the student, plus each 
student's share of the variable costs attributable to faculty. Tui- 


tion for two years and a summer comes to $6750; allowing $200 a year 
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and 350 a summer for chemicals and laboratory specimens, computer 
budgets, and the like, leaves $6300. This is the student related 
portion of To, Tz, and Tq. 

Primary teaching teams are presumed to be comparable in cost 
from one program to the next, though their professional composition 
will differ somewhat. The representative team for the Health Asso- 
ciate program includes a physician, at $27,000 annually, two social 
generalists at $15,000 annually each, and three health practitioners, 
at $13,300 each annually. Each such primary teaching team is sup- 
ported by 1.3 secretaries, at $6800 per secretary per year. This 
comes to $105,740 in direct salary expenses. Adding 14%, the cur- 
rent Johns Hopkins fringe benefit figure, plus $200 per secretary 
per month for office supplies, plus $12,000 per team for media 
resource support, yields $135,664, the total variable cost per 
teaching team. Then, over the two year curriculum, the Health 
Assistant student is responsible for 2/30 of the annual costs of one 
such team, or $9044, and each student of the other programs is 
responsible for 2/25 of these costs, or $10,853. 

Faculty specialists do not contribute to the T; because--as 
noted earlier--in the judgment of the staff of the School, their 
numbers will not vary significantly. 

The final contribution to the T; is the per student costs of 
additional clinical supervisors. The numbers of days weekly that 
students in each year of each curriculum will receive clinical 


instruction is shown in the following table: 












Prae a e ሙሄበሪመጦ 2 


መ መጨ w መ  ፒ ጡወ‏ سي ې 





<A -—— -— o — له‎ ዘብ መሙ፡-*›-8 
— س — سم وھ تا‎ -—— 4, 0 
— c سم ومس‎ A 0 A -- o» am ا‎ 


899(ሸበ9፡6፡፡5569966695 A =  ‹ኡ፡ e e 







om weer (ee ص‎ መሙ 
DE NE ہب نګ‎ o M 1 OM سیف‎ 
08 ሙተ HS > 4  — UU. ٢ 


| — o—e — o — eee لن و‎ - 


6 6 . سا‎ ሙ - "E 
— Geo cnn em 
E  —— ٩1 بپ‎ a P 
= t- [> مس‎ LI 3 سے‎ 
۶ص — سص أ‎ A — ' 
= ٩ ጠው 0 © > መሙ - Y 
a ጠሬ ces سم‎ Å am سے‎ 

E (65 በጠ - "- we p‏ بب 


68 


H. Ass't H. Ao/N.P. H. Serv. Mgmt. Env. Hygiene 


program prozrams program program 
18% year 
students ፲ 5 I hour 2 hours 
2nd «ear 
students 35325 25 1 hour 2 hours 
Sum 4 6 ለ ን) 0533 1068 


Clinical supervision will be carried out in groups containing ons 
instructor and two students, and each primary teaching team will 
perform the needed supervision of 8 of its own students (two for 
each of the four members qualified to perform clinical supervision). 
Then the average requirement (in days weekly) for additional clini- 
cal supervisors is 22/30 of the Health Assistant program figure 


above, or 17/25 of the figures for the other programs: 


Health Assistant program: 5117 
Eealth Associate and Nurse 
Practitioner programs: 2.94 
Health Services Management 
programs: 40362 
Environmental Hygiene 


programs: 0725 


Each clinical supervisor instructs two students at once, 5 days 


weekly; then his capacity is 10 student-days per week. Dividing the 
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previous row by 10 gives the following numbers of additional clini- 
cal supervisors required per student matriculated into the indicated 


curriculum: 


Health Assistant programs: 7 
Health Associate and Nurse 
Practitioner programs: 0294 
Health Services Management 
program: «0036 


Environmental Hygiene program: 2 


Full time equivalent clinical supervisors will be either physicians, 
salaried at roughly $27,000 annually, or health practitioners at 
about $13,300 annually. On the assumption that their relative pro- 
portions will be roughly 18 of the former to 10 of the latter, the 
average direct salary becomes $22,107; adding 14% for fringe bens- 
fits yields a net cost of $25,202 per additional clinical super- 
visor. It is presumed that the office and secretarial support 
requirements of additional clinical supervisors are negligible. 


Vultiplying the figures in the previous line by this sum gives: 


Health Assistant program: $7855 
Health Associate and Nurse 


Practitioner programs: $7409 
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Sumarizing, the income less proportional costs of instruc- 


tion per matriculated student are shown in the following table: 


Health Assistant H. Associate/ H. Serv. 


program N. Practi- 
tioner 
progrems 
Direct student 
revenues less 
direct student 
costs $16,899* $6300 
Primary teach- 
ing team costs $-9044 $-10853 
Additional 
clinical super- 
visor costs 53-7855 $-7409 
Total $0 $-11,962 


*Includes $10,599 from the NIH. 


۲ ۲ 0۵ 2 


$6300 


$-10853 


$-91 


$-4644 


Env. 
Hygiene 
program 


36300 


$-10853 


$-181 


$-4734 
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IV-4 CONSTRAINTS 

Two of the constraint families listed and discussed in 
section III-3 are encorporated into the School of Health Services 
equilibrium planning model; they are faculty availability and 
facilities availability. The other listed constraint types are 
inappropriate either because they are obviated by the construction 
of the model (e.g. the mathematical identities relating faculty and 
students, or matriculations and enrollments), because there is no 
foreseeable prospect of violating the availability of the resource 
in question (student applications), or because the constraint is 
simply not an issue for the School (faculty or student body composi- 
Cion). 

The only faculty category which in the long run is antici- 
pated to be in short supply is part time clinical supervisors. 
These will be called "additional clinical supervisors", for their 
supervision is in addition to that provided by members of the pri- 
mary teaching teams. Thus the model includes one faculty availa- 
bility constraint. The facilities availability family, on the other 
hand, contains eight constraints because eight different facility 
types which might conceivably be in short supply have been identi- 
fied. These are faculty offices, three types of laboratories, and 


four types of classrooms. 


Availability of additional clinical supervisors: 


It appears trat the only faculty category whose availability 
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might possibly be limiting is additional clinical supervisors. The 


appropriate constraint is 


where aq]; is the number of additional clinical supervisors required 


per student matriculated in the iP 


curriculum and A, is the number 
of additional clinical supervisors realistically available to the 


School, in full time equivalents. The estimates for the a); are: 


11 = .4033 
8712 5 8 

213 = 212130052 
84 = 5 


These were developed as follows. It is presumed that clinical 
supervision is carried out in groups of one clinical supervisor and 
two students. Of the six member primary teaching team, four members 
will provide clinical supervision for students of their sections; 
since &ll members of a section undergo clinical supervision at the 
same times, this means that the primary teaching team can accommodate 
the clinical supervisory requirements of eight students, leaving 22 
students from each health assistant section or 17 students from each 
section in other curricula who will need to be accommodated by addi- 


tional clinical supervisors. Clinical supervision is most intense 
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during the spring terms; thus it is on the requirements at those 
times that the estimates are based. The expected spring term inten- 
o? clinical supervision in days weekly for students of each 


year of each curriculum is shown in the following table. 


H. Asstt. H. Ao/N.P. H. Serv. Mgmt. Env. Hygiene 


18% spring term 1 1.5 1 hour 2 hours 
2nd spring term 5 4.9 1 hour 2 mours 
Sun 5.9 6 00533 » 1066 


Multiplying by 22/30 or 17/25 as appropriate yields the average num- 
ber of days weekly each such student requires an additional clinical 


Supervisor. 


H. Ass't. H. Ao/N.P. H. Serv. Mgmt. Env. Hygiene 


4.035 4.08 00362 0725 


Bach full time equivalent additional clinical supervisor instructs 
two students at once, five days weekly; his capacity is 10 student 
days weekly. Then dividing the previous row by 10 gives the number 
of additional clinical supervisors required per graduate of each 
program; these are the figures cited earlier. 

A,, the number of additional clinical supervisors available, 
in full time equivalents, is estimated to be 28. This is based on 


the presumptions that each interested individual contributes 1/10 
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full time equivalent--i.e. a morning or an afternoon a week, that 
25% of the Johns Hopkins house staff in Pediatrics, Medicine, 
Surgery, Psychiatry, and Obstetrics (totalling 399) will be inter- 
ested, and that 10% of the general practitioners (420), internists 
(654), vediatricians (399), and OB-GYN specialists (365) in the 
Baltimore metropolitan area will be interested. The hospital staff 
figures are as of September 1972; the Baltimore metropolitan area 
figures are from an American Medical Association survay report for 
the end of 1969, the last year in which the AMA cumulated and pub- 
lished its survey by metropolitan areas. Apparentiy hospital house 


steff are not included in the survey. 


Availability of faculty offices: 


The appropriate constraint is 


Pm» 


5, 82; 04*21, = Ay 

1-1 
where as; is the number of faculty desks required ver student 
matriculated in curriculum i and where An is the total number of 
desks available. The possibility exists that offices might need to 
be converted to dry (instrumentation) laboratories; the term 2Y} 
reflects this drain upon the number of available desks, assuming 
there are two desks ver office. It is presumed that only the vri- 
mary teaching teams and additional clinical supervisors are provor- 


Gional to enrollments, and that of these, only primary teaching team 
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members will require desks. In the Health Assistant program, each 
student's share of primary teaching team members is 1/30 ور سیر‎ 
each year, or .4 for the two years; in the other curricula, it is 


1/25 of 6 for each of two years, or .48. Then 


821 = 4 


A‏ سس پس د 


Hampton House, formerly a student nurse's residence, is cur- 
rently undergoing a prolonged renovation; when complete, this build- 
ing will house the School of Health Services and portions of other 
divisions of the University. Renovation biueprints for the floors to 
be occupied by the School show space for 80 desks for teaching 
faculty. Approximately ten of these will be required for selected 


faculty specialists; then An = 70. 


Facilities usage informations: 

Six types of instructional facilities (three sizes of class- 
room, and three types of laboratory) have been identified as poten- 
tially being in short supply. The classroom types are large class- 
rooms--that is, classrooms capable of holding a double section, con- 
sistinz of 50 to 60 students, depending on the curriculum; medium 
classroons--i.e., classrooms capable of accommodatinz a single 
section but not two; and small classrooms--i.e., those capable of 


holding a half section but not an entire section. Bench laboratories 
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are "wet" laboratories fitted with laboratory stations for each per- 
son; these labs are of the sort used for practical instruction in 
chemistry, biology, etc. It is anticipated that entire student 
sections will perform their bench laboratory studies at once; thus 
the capacity of a bench laboratory must be 25 to 50 students. 
Demonstration laboratories are wet laboratories in which the student 
will observe laboratory experiments and procedures performed by an 
instructor. Demonstration laboratory sessicns will be attended by 
half sections; therefore, the capacity of each demonstration labora- 
tory should be at least 13 to 15 students. Dry laboratories are 
rooms furnished with special equipment required in practical 
instruction, for example a bed on which one student may recline 
while another takes his pulse and blood pressure. It is anticirated 
that dry laboratories will be utilized by groups of four students. 
The following table contains the numbers of hours weekly 
that students of each curriculum are expected to spend under 
instruction in the indicated type of facility. The figures are for 
the fall semester of the year because usage of these facilities is 
heavier during the fall than during the spring. The symbol I 
denotes the fall of the first year of the indicated curriculum; II 


denotes the fall of the second year. 





Hi Aes't. 


H. Assoc/ 
KOTEREC. 


H. Serv. 
M gmt ወ 


Environ. 
Hygiene 


۱ 


II 


1890 


1۳1 


Larg 
2 


2 


Classrooms: 


e Medium Small 
5 5 
5 4 
x 5 
3 5 
225 3-15 
25 375 
15 ۱۱ 20 
7 3.75 3 ۶5 


1 


Laboratories: 
Dry Demo. Berch 
1 2 1 
1 2 1 
Delo 92 0 
0 1 1 
3:213 0 0 
379 0 0 
1.87 1.97 ETS 
0 0 0 


FACILITIES USAGE TABLE 


(hours weekly) 
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Dry (instrumentation) laboratory availability: 


The appropriate constraint is; 


y e3103-37.5Y] É Ag 

0 
where ag; represents the number of hours (weekly) of dry laboratory 
availability required per student matriculated in tne jth curriculun, 
and Ag represents the number of hours of dry laboratory usage avail- 
able weekly. The term -57.5Y,, which is based on a 7.5 hour school 
day, represents the contribution in hours per week of each rocm 
which is nominally a faculty office but which is converted to dry 
laboratory use. 

Dry laboratory useage rates per student are the figures from 

the facilities usage table summed over the length of each curriculum 
and divided by four (reflecting the expectation that dry labora- 


tories will be used by four students at once). Thus 


8 6 
Azo ኖ «95 
855 # 1.875 
arg = ,468 


The availability of each type of laboratory is not presently 
mown. No laboratories of any kind are planned for Hampton House; 


excess capacity elsewhere on the Johns Hopkins East Baltimore campus 
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is difficult to assess. It is not generally possible to estimate 
excess laboratory capacity by tallying up unused rooms, because 
excess capacity usually manifests itself in a spreading out of 
ongoing activities and in creation of new ones which would not have 
been undertaken at all had "free space” not been available. Never- 
theless, it is known that some facilities will be made available to 
the School, and certainly it can be assumed that the closure of ths 
Hospital's nursing program will result in some free space. A 
reasonable minimum estimate is one laboratory of each type continu- 


ously available. 
Then Az = 37.5. 


Demonstration lacoratory availability: 


the appropriate constraint is: 


4 
? Le 
Di. دا زره‎ * A4 

i=] 


where aq; is the number of hours weekly of demonstration laboratory 
usage availability required per student matriculated in curriculum i, 
and where A, is the total availability of demonstration laboratories, 
measured in hours weekly. 

Demonstration laboratory requirements per student are the 
figure from the facilities usage table summed over the length of the 
curricula and divided by 12.5 or 15, reflecting the expectatation 


that demonstration laboratories will be used by half sections. Then 
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a4] = 6 
AE 
843 = O 

24 = «150 


The difficulty in estimating available laboratory space has 
been mentioned in conjunction with the previous constraint. It is 
felt that at least one full time demonstration laboratory or its 
equivalent in part time availability will be available; then 


Aj = 3745. 


Bench laboratory availability: 


The appropriate constraint is 


Tm dl 
و‎ 6+۳۷ 
1-1 


where asi is the amount of bench laboratory usage required (in hours 
per week) per graduate of curriculum i, and where As is the expected 
availability of bench laboratories, in hours per week. 

Bench laboratory requirements per student are found by sum- 
ming the figures in the facilities usage table over the lengtn of 
each curriculum and dividing by 25 or 30 as appropriate (recall that 
bench laboratories will be occupied by entire student sections). 


Then 
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as] = 4.0867 
Sco, "UMS 
ape = ۵ 

Sca. 7 و‎ 


Again, in the absence of reliable information on theswiti- 
mate availability of laboratory space, one full time laboratory is 


presumed; As = 37.5. 


Large classroom availability: 


The appropriate constraint is 


4 


u 25:4 *Y2*Y5 * Ag 
where ag; is the requirement in hours per week for large classroom 
space per student entering the ¡th curriculum, and where Ag is the 
availability of large classrooms in hours per week. The terms +Y, 
and +Y reflect the effects of using large classrooms for classes 
other than (combined) double section classes (recall that Yo is the 
number of hours weekly that large classrooms are partitioned and 
used for smaller classes, and Ys is the number of hours weekly that 
large classrooms are used for single smaller classes). 

large classroom requirements are found by summing the figures 
from tne facilities usage table over the length of each curriculum 


and dividing by 50 or 60 as appropriate; then 
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op ከዱ. 
ago = 68ء‎ 
agz = 550 
aga = 7ھ‎ 


Hampton House renovation blueprints show one large seminar of 
approximately 630 net square feet (this room is to be equipped with 
a divider) and three of 560 net square feet each; any of these four 
classrooms could accommodate a double section (though perhaps not 


very comfortably). Then Ag = 4 x 37.5 = 150 hours weekly. 


Medium classroom availability: 
The appropriate constraint is 


۶ سه 
1 


ረ 
۱ M Qu 


4 
01-7 


where aq; is the requirement in hours per week for medium classroom 


Secs per student entering the ¿th 


curriculum, and where Av is the 
availability of medium classrooms in hours per week. The term - ولغ‎ 
reflects the effect of using a partitioned large classroom for two 
medium classes simultaneously; the term - Yz reflects using a large 
classroom for a single medium class. The last term, *Y,, incorpor- 
ates the utilization of mediun classrooms for small (half section) 


classes. 


Medium classroom requirements per student are found by 
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suminz the figures in the facilities usage table over the curricula 


and dividing by 30 or 25 as appropriate. Then 


a71 = 3 
272 = 4ه‎ 
azo = 1.8 


Hampton House renovation blueprints show one seminar of 
approximately 400 net square feet, two of 390 net square feet, 
another of 360, one of 330, and two of 315 net square feet each. 
Each of these seven rooms is capable of accommodating a full section; 


then Ay = 7 x 57.5 = 262.5 hours weekly. 


Small classroom availability: 


The appropriate constraint is 


ya 
agiCi-Ya ” Ag 


ولا 
3 .1 


where ag, is the requirement in hours per week of small classroom 


availability per student entering the jth 


curriculum, and Ag is the 
availability of small classrooms in hours per week. The term -Y, 
reflects the relaxation of this constraint which occurs when Y, 


weekly hours of medium classroom availability is devoted to small 


classes. 
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omall classroom requirements per student matriculated in 
program i are found by summing the figures in the facilities usage 
table over the length of each curriculum and dividing by 12.5 or 15 


as appropriate. Then 


agi x vb 
ago x .8 
8 ae <0 
aga 1٤8 


Hampton House renovation blueprints show one seminar of 
approximately 210 net square feet; this room will accommodate one 


half student section. Then Ag = 97.5. 


Availability of partitionable larze classrooms: 

During the discussion of the large classrooms availability 
constraint it was noted that one of the large classrooms is to be 
equipped with a divider; then the amount of partitioned large class- 
room usage, represented by the variable Yo, can not exceed the num- 
ber of hours available weekly in this one partitionable room. Then 


the appropriate constraint is 


۳4 
Yo ۹ $7.49 





IV-5 THE MATHEMATICAL MODEL 


The following is & mathematical summary of the School of 
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Health Services equilibrium planning model, preserted, for conveni- 


ence, in matrix format: 


Subject 


€ 


ወቁ: 


| 
| 
۱ 
| 
| 


Maximize 

to 

. 408 00362 
8 8ء‎ 
.94 1.875 
228 
DIS oO 
028 0 
24 8 
.8 .6 

0 0 


1 0 0 0 

0 1 0 0 

0 0 1 0 

0 0 0 1 

O -11962 -4644 -473£ 
0725 0 0 ل‎ 
8 2 0 0 
98 -37.5 O 0 
.150 0 O 0 
.15 0 OTTO 
mw o O 1 1 
075 0 -2 => 

ار 0 2 1 
0 1 0 0 
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The results of the actual manipulation of this model will be 


presented in Chapter VII following a detailed presentation, in Chap- 


ters V and VI, of the general method of finding the efficient mani- 


old 
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NETHOD OF SOLUTION 


V-1 INTRODUCTION 
This chapter is concerned with finding solutions to the 
following general problem, which will be called "problem P"; 
Find the set of efficient solutions to 
n 
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in 
“21 = 0 
“22 = 0 
Su res 0 


where it is assumed that the set of all feasible solutions to 
the constraint set is bounded, and that there do not exist 


positive scalars L; such that 


3 
2, 54541 " 0 for j 7 1,2,...n 


It was noted in Chapter III that efficient solutions to this 
problem are of interest because the preferred solution must be among 
them. In this chavter, a technique for finding the complete set of 
efficient solutions to problem P will be presented. The reader with 
no taste for algorithms or mathematical proofs is invited to skip to 
Chapter VII where he will find an account of the application of this 
technique to the School for Health Services equilibrium planning 


model wnich was developed in Chapter IV. 
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one procedure for rinding the efficient poinmvsewill ve» pre= 
sented in two parts. First, an algorithm for finding all efficient 
vertices of the polyhedron of feasible solutions will be presented. 
It will be demonstrated that the set of all efficient solutions can 
be expressed in terms of these efficient vertices. Then the second 
part of the procedure, an algorithm for finding such an expression, 


will be presented. 
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V-2 BASIC DEFINITIONS AND NOTATION 
The following notation and basic definitions will be em- 


ployed in the remainder of this chapter and in Chapter VI. 


Inequality notation: Let x = رو ہے وی ہیا‎ and y = وت‎ 


be two vectors in n-space. 


Cee x 


(2) x 


ly 


7 x; = y; for Sul 1 = lu 1 


ir 


Voit x = y and x; > y; for at least one i. 


SN ۳955۱۱۲۷۰۰۱۰۶7 


Transpose notation: When t is employed as a superscript it will 


indicate the transpose operation on the superscripted entity. 


Inner product notation: Without ambiguity, xy will signify simple 
multiplication - i.e. "x times y" when x or y or both are 
scalars, and will signify the inner product of x and y when 


both are vectors or matrices. 


Components of vectors: Components of vectors (and rows or columns 
of matrices, when specified) will customarily be represented 
by subscripting the symbol for the vector (matrix); e.g., x, 
is the it component of vector x. Occasionally, components 
of the results of? inner product operaiions will be referred 
to in the same fashion. For example, if L is a vector and G 
15 a matrix of appropriate dimension then the j 6ت‎ 0980-0 


of the vector resulting fron the inner product of L and G may 


ba written (LG);. 
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"1" and "0": It is sometimes necessary to employ a vector every 
component of which is one, or a vector every component of 
which is zero. The symbols 1 and O will be used, respec- 
tively. In every case the information in the remainder of 
the equation or expression will be sufficient to infer 
wnether these vectors or the scalars zero and one are 
intended, and in the former case the length of the vectors 


will be clear also. 


Problem P (vector and matrix notation): 


Define: x = the column vector of variables x; . x has 
length ne 
L= the row vector (of length k) of multipliers Li.» 
C = the matrix of criteria coefficients C has 


“41 
k rows and n columns. Notice that the coeffi- 
cients of each particular criterion appear as 
a row of C. It will be assumed that there 
does not ezist L > 0 suchthatile = 0; 

F = tha "feasible region" - that is, the set of all 


x satisfying the given constraints: 


n 
Jl 
P ረ 
: e e e - . 1 = + © a ቀ 
0 RD ¡opt 1 ووو‎ 
7 ۲ 
ን د‎ ISO] ANM 


t is assumed that F is bounded. 
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the matrix of coefficients on the left of tne 


='s in the equation set above. A has dimen- 


sions (mtn) by n. 


the column vector (of length mtn) on the right 


of the ='s above. 


In this notation P is simply: 


Find the efficient points to 
(Max) Cx 


subject to x in F 


Feasible solution x will be said to be efficient to 


problem P if and only if there does not exist y in F such 


that 


iv 


The terms "vertex", "extreme point”, and "corner" will be 


used interchangeably to denote the solution to a "regular 


system of tight constraints" (see any linear programming 


text). 


Remark: 


It is well known that any point in a bounded 
polyhedron (such as F) may be expressed as a 
convex combination of its vertices. For 
example, if x is in F and if xt,x*,...,x¥ are 
thə vertices of F, then there exist multipliərs 


ns such that 
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Also, the intersection of a polyhedron and of 
one of its constraints taken as an equality 
constraint is another polyhedron, whose ver- 
tices are vertices of the original polyhedron. 
(The geometry of polyhedra is well explained 


in Simmonard.!) 


~ 


Tight constraint: A constraint (whether originally stated as an 
equality or as an inequality) which is satisfied as an 


equality is said to be “tight”. 


Facet: A facet is the intersection of a polyhedron and one or more 
of its tizht constraints. 
Remark: It may be trat a facet is vacuous = if, for 
example, the chosen tight constraints do not 
intersect in F. By the remarks above, vhen a 
facet is nonvacuous, it is a polyhedron every 
point of whicn may be expressed as a convex 
combination of its vertices - that is, of some 
of the vertices of F. Clearly every nonvacu- 


ous facet arising frm F must contain at least 
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one vertex. The symbol f(J) will denote the 
facet consisting of the intersection of F and 
the tight constraints whose indices are in the 
set J. Notice that if set J contains set ال‎ 
then facet f(J') contains facet f(J), for f(J) 
is the intersection of f(J') and those tight 
constraints indexed in J but not in J'. Inci- 
dently, vertices are facets, for they are 
intersectionsof F and of a regular system of 


tight constraints. 


Efficient facet: An efficient facet is a facet of which every solu- 


tion x in the facet is efficient. 


Efficient manifold: The efficient manifold is the set of all effi- 


cient solutions to problem P. 
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V-3 PRELIMINARY THEOREMS 

The literature of multicriterion optimization was reviewed 
in Chapter I: the fundamental results on which the efficient point 
search algorithms are based are presented in this section. 

Problem P includes the assumption that there does not exist 
L > O such that LC = O., The following theorem establishes that the 


provlem is uninteresting when this condition is not met. 


Theorem V-1*: Every x in F of problem P is efficient if there 


exists L>" Oi suchi that ICE = O, 


2 


Proof: By Steimke’s theorem of the alternative” the existence of 


row vector L > O such that L = O implies that there does not 


ከ 


exist x such that Cx £ O. Since LC = 0, then L(-C) = O, and 


لاا 


there does not exist x such that Cx > O. Consider any two 
solutions x' and x" such that x’ # x". Define x = x"=x'. 
It has been seen that Cx jf O and Cx f O. Notice that Cx £ O 
implies either 

a) Cx = Û, or 

b) (Cx); > 0 for 660 
Also, Cx 7 O implies either 

c) Cx = O, or 


d) (Cx) j < O for some je 


at. ; : 

Throughout this chapter and the next, theorems, lemmas, and 
corollaries not expressly credited to someone else were derived by 
the author and are, to the best of his knowledge, original. 
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Together these imply either 
e) Cx = 0, or 
f) (Cx); د‎ 0 for some i and (cx); < 0 for some j. 
That is (fram x = x"-x') either 
ey Cx Cx’, or : 
h) (Cx"); > (Cx'), for some i and (Cx") 5 < (Cx') 5 for 
soma j. 


Then Cx" P Cx'. Since x" was arbitrary, x* must be effi- 


ማከር 


cient. 


Since x’! was arbitrary, all of F is efficient,» 


QED 


It happens that the existence of L > O such that LC = 0 


iV 9 


implies and is implied by the existence of L = 1 such that LC = 0. 
Let vector M= L - 1 and nove that the assumption may be reexpressed 
"there does not exist M = O such that MC = -LC". This condition is 
in a form easily tested by phase one of the two phase simplex 
method. 

The following two theorems date in principle back to 1951, 
to Kuhn and Tucker's fundamental paper on nonlinear programming. 
These, and corollary V-3-1, are cited and proven by Philip.? Problem 

"ncn 


P' is defined to be problem P modified to include only “= con- 


straints. 


Theorem V-2 (Philip): Let w = LC for some L > 0. If x solves Max wx 


subject to x in F, then x is efficient to problem P'. 


| nn 
6 - 
۳ 


AA fi P9 ae—. P اح‎ 


N Eu ے - وت اد‎ 
L1.  ——— 11 1 IDE l 


09 - ከ. مھ‎ 
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Theorem V-3 (Philip): Let J be the set of indices of the rows (a; ) 
of matrix A of problem P' which are tight at some solution 
Rew Lite x TS Er ICIC natto proulem F'’™ar and only if there 

> 


exist row vector L > O and multipliers عا‎ = O (for all i in 


J) such that 


=e = pi a; kK; 
ali IRS 
Colollary V-3-1 (Philip): Under the conditions above, x is efficient 
to problem P' if and only if there exists L = E > 0 for 


arbitrary vector Eş and k; = 0 for all i in J such that 


I 
ا‎ > ), 29.51 
a rr) 


The following ccrollary extends this necessary and sufficient condi- 


tion to problem P. 


Corollary V-3-2: Let J be the set of indices of tight constraints 
at x in F of problem P. Solution x is efficient to problem 
P if and only if there exist L = E > 0 and k: for all imn 
J, where k; = O if constraint i is an inequality and unre- 
stricted if constraint i is an equality in the statemant of 


F, such that 


LC = > asks 


quy 
all i in J 





Proof: 
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(If clause) The equations defining F may be grouped into an 
equality set and an inequality set; letting Al pt ae, and b? 
be tha coefficient and right hand side matrices and vectors 


of these respectively, F may be expressed: 


Aix = pt 


Ax = be 


An equivalent expression for F is: 


ኤ 
= 
ዕና 
11 
c 
مم‎ 


>‘ 
دم 

۷ 
ዘኑ 

er 
دم‎ 


ዘኑ 


a 


With F specified in this all inequality manner, the problem 
is ready for application of corollary V-3-1l; that is, solu- 
tion x is efficient if and only if there exists row vector 

LSE مه‎ nonnegative multipliers k'; for all i in J(and 


in the third set of equations above), such that 


LC = /, ak’. 


Qe 
all iin J 


These k,' may correspond to equations in any of the three 
groups of equations in the inequality expression for F. Let 
J be partitioned into sets 31,38, and JÍ accordingly. Then 


the previous expression becomes: 


وو مه 
p መሸማ‏ 


+ 
= ا هه 


nd. پسم‌ی.‎ es وټ‎ 
De ئ908۵ ه َ يه مي‎ 
_—_ «2 መ- o o -——— c8 ‹> 
SS ARA *ዬዒያ 

5c c—— 6 
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isind iin له‎ 96 118 0 


All equations in the first and third groups of equations in 
the inequality expression for F must necessarily be satis- 
fied as equalities, so for every a; in the first group, 
there is a term in the first summation of the last equation, 


and also another for موہ‎ 


1 in the third summation. Grouping 


these two, making the sign substitution, and representing 


ki MEO J$ ከሃ ki" gives: 
وت‎ 
LC= رگ‎ .s4(ki'-ky")* Zs او‎ 
iin J i in J* 


Setting k, = (k,'-k,") from the first group and k, = k,' for 


J? gives a set of k* satisfying 


LC = ረ az ks 
Jom: d 
= ; ደ : : 
ks = O when constraint i is an inequality 


This is the condition specified in the statement of the cor- 
ollary. 

(Only if clause) Since every real multiplier k; can be 
expressed as the difference of two nonnegative multipliers 
k,' and k,", a set of kj satisfying the condition in the 
statement of the theorem can be used to create a set of k;' 
satisfying tne corresponding conditions for the all inequali- 


ty version of prooulemer’. 


QED 
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The following definitions and concepts are familiar from 
linear programming; the only novelty is the occasional addition of 
an extra subscript reflecting the fact that the cost vector is now a 
cost matrix. Let: 

B = any basis (later on, bases will be identified by sub- 

scripts, e.g. Bi) 


سل 
c; = the jh row of objectives matrix ©‏ 
the vector of those components of c; corresponding to‏ = وله 


the basic components of x for the basis B. 


5.47” the jth column of matrix A 
- pol 
Me = B دہ‎ 
Y = the matrix of columns Y3 
2 B 
> = AC iJ j 


s ከ This is the reduced cost of the ¿Eh 


component 
of x calculated with the current basis and the i"! pow 
of C, 

G = the matrix of reduced costs Sij 


- 


= 8”? This is the vector of values of the basic’ varia- 


a 
t 


bles. 
The following tableau, the usual simplex tableau with a 
reduced costs matrix instead of a reduced costs row, will be 0 


the "main tableau"; 
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5 5 
yo.‏ د ان ا 


Se ገ] 


۷ | 


Figure III 


Philip mentions without proof that the existence of L > O 


ዘሎ 


such that LG = 0 is an efficiency criterion, analagous to the opti- 
mality criterion of linear programming. This condition will be 
referred to as the “efficiency criterion”. Lemmas V-4 and V-5 and 


theorem V-6 establisn the validity of this criterion, and expose an 


idiocyncracy unmentioned by Philip. 


Leama V-4: If there exists row vector L > O such that LG $ O for 
the matrix G corresponding to a tableau representing solu- 
tion x to problem P, then x is efficient to problem P. 


Proof:  Follows directly from theorem V-2. 


Lemma V-5: If x is efficient to problem P then there exists vector 
L >0O such that x solves Max LCx subject to x in F. 

Proof: By the corollary to Philip's third theorem, x is efficient 
implies there exists L > 0 and multipliers k, ( = 0 when 
constraint i is an inequality in the statement of problem P) 
such that 

ہے 


= 3 2 ۰ 
LG zer: Ay Ky 
1T ud 
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where J is the set of constraints which are tight at x. The 
vector L satisfies the requirements of the lemma, for if it 
did not, then there would exist x' in F such that LCx' > LCx. 
Substituting the relation 
); 
LC = ask: 
"0 
into the above expression yields 


د در سج ee‏ 


1 in J a in J i in J 


en since ajx — bi fem i ind, since Y is the set of tight 


constraints at x, 


= سے 
0ج ንን k; (83x')- s ks (b4)‏ 
18 1 مد بد 


or 


De k,(a;x'-b,) 9 10) 


i in J 
Kk, may be negative only when constraint i is an equality; 
then a,x'-b; = 0, so no term containing k; £ O contributes 
to the sum abovs, Then there must bs at least one term with 
ks » O and a,x'-b, > O which is a contradiction since x' is 


feasible. 


QED 


V-6: If solution x is an efficient vertex to problem P, 


there exists a basis B corresponding to x such that for the 
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associated matrix G there exists vector L > O such that 
jos 

Proof: By the previous lemma, there exists vector L > O such that x 
solves max LCx subject to x in F. Ina concluding tableau of 
this problem corresponding to x the reduced costs row is 


nonpositive. That is to say: 


or 


ቸን ን ከራ 
کے کا‎ )y; = Oven ALIN 
qe 


C 


Tg 


i3 
i=1 J 


or 


m 5 1 
2۳ 1 y 5) = 0 forall j 
17 


or, using the relation a کر‎ 54.74 
= Z 
T = 0 for all j 
ግ በከ ل‎ 
or 
LG ٤ ٥ 


QED 


These results are not the same as saying that x is efficient 
if and only if such L exists, for when the vertex at x is degenerate, 
there are several corresponding bases and corresponding reduced 
costs matrices. The last theorem proves only that one of them must 
satisfy the efficiency criterion; an example containing s degener- 


ate vertex for one basis for which such an L > O does not exist is 
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given in Appendix I. Then the existence of such an L for a given G 
is sufficient but not necessary for the efficiency of the vertex 


represented. 
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V-4 AN ALGORITHM FOR GENERATING ALL EFFICIENT EXTREME POINTS 

Theorem V-2 provides an easy method for finding an initial 
efficient vertex. The procedure is simply to take some arbitrary 
strictly positive weighting of the rows of C, using it as the objec- 
tive gradient of an ordinary linear program (to be maximized over 
F). The simplex will terminate at a vertex, which must be efficient, 
according to the theorem. 

Philip® presents a method, based on the simplex tableau, for 
determining whether a given efficient vertex is the only efficient 
vertex, and if not, for determining anotner vertex, one which may bs 
reached by one pivot. (A vertex or basis which may be reached from 
another by one simplex pivot will be said to "neighbor" the first). 
۲ ٥ 6> alcorithw! to follow, this procedure will be used to find all 
neighboring efficient vertices of the initial vertex, and then all 
neignbors of those neighbors, and so on. It will be proven that all 
efficients vertices can be generated by this procedure. 

Suppose that the linear program "Max LCx subject to x in F” 
is solved for some fixed positive vector L. At termination there 
will be a basis B for which the associated reduced costs matrix G 
has the property that there exists L > O such that LG £0. It is 
desired to know whether any neighboring vertices are efficient also. 
The following procedure for answering this question is due to 
LID 

suppose that there exists vector L' such that L' > 0, 


ፌ 
EG 0 and (L'G); = Q0 for some i. Then both the vertex represented 
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by the present tableau and that represented by the tableau which 
 —— if column i were pivoted into the basis solve the prob- 
lem "Max L'Cx subject to x in F'". (Clearly the current basis solves 
this problem, for L'G É O implies the linear programming optimality 
criterion is met; when column i is pivoted in, since (L'G), = 0, the 
linear programming reduced costs row is unchanged; the new tableau 
meets the optimality criterion also). 

It happens that the new tableau may be simply another 
expression for the same vertex. In this case, the vertex is degen- 
erate, and it can happen that not all neighbors of a degenerate ver- 
tex are uncovered in the course of investigating all bases neighbor- 
ing a given basis. Nevertheless, the procedure outlined will in the 
end find all efficient vertices, as theorem V-10 will show. 

Thanks to the homogeneity of the requirements L > O and 
2 O, if such L exists then there exists L=E>0 for arbitrary 
positive E, and vice versa. Then an equivalent efficiency condition 
is the existence of L = 1 such that LG £0. Let HM. ማም 8... 
perform a change of variables; then the condition becomes: the ver- 
tex represented is efficient if there exists M= 0 such that 
MG $ 1G. This condition is easily checked through phase one of the 
two phase simplex. In checking the efficiency of a neighbor, it is 
desired to find such an M with the additional property that 
(MG), = (16); for the particular column i in question. This, also, 
is easily checked; simply minimize the value of the slack to the 


constraint (NG); = (1G); in the problem just outlined; if the problem 
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is feasible and the slack can be driven to zero, then the desired 
condition is met; the basis resulting from the introduction of 
column i will satisfy the efficiency criterion. 

Since it is desired to find all efficient neighbors, several 
different pivot locations may need to be investigated for any given 
column i satisfying the entry criterion just discussed. The exit 


criteria of interest are the following: 


(a) وذ لاد‎ - minimum (a. /y. i) 
s/Y م تې‎ 0 


(d) d, = O anà Yj FO 


The reader may verify that for column i satisfying the entry cri- 
terion just discussed, any pivot row satisfying either of tne condi- 
tions (a) or (bd) above will give a new tableau whose reduced costs 
row (for the objective L'Cx) will satisfy the linear programming 
optimality criterion--thus the associated G satisfies the efficiency 
BriLtericn. 

The following lemmas and theorem V-10 establish that a 
neighbor by neighbor search of the type suggested will generate all 


vertices efficient to problem P. 


Lemma V-7 ("Connectivity Lemma"): Let Bj and Bo be two optimal 
bases solving, for fixed L > O, tne proolem Vax LCx subject 
to xin F. There exisus a Pate sequence ias ehboring 


bases--whose tableau meet the efficiency criterion--linking 
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the two. 

Proof: It is well known that the simplex may be used to find all 
alternate optima to a linear program. In so doing, it 
passes from basis to neighboring basis (by pivoting accord- 
ing to the exit criteria (a) and (b) above), each of which 
satisfies the linear programming optimality criterion--and 
hence the efficiency criterion. Suppose that the process 
begins at By; eventually Bo is found. Tracing backward 
along the path of discovery yield the sequence sought. 


QED 


The following rather pedestrian identity is required in 
later proofs. Parametric objective functions will figure in these 
proofs, but for the moment it is sufficient to understand that V' 
and V" are two objective gradients (that is, V'x and V"x are objec- 
tive functions) and a third, V(s), is given by V(s) = V' + s(V"-V') 


where s is a scalar. 


Lemme V-8 ("Reduced Costs Identity"): Let V' and V" be two objective 
gradients and let V(s) » V' * s(V"-V'*). Let E19 be the 


reduced cost for the variable x, corresponding to objective 


J 
V(s) and basis Bi, and let g(j,i) be the reduced cost for 
that variable corresponding to basis B, and the objective 
(v"-V')x. Then for any given scalar s*, 

e(j,i,s) = g(j,i,s*) * (s-s*)g(3,i) 


Proof: By the definition of reduced cost, 
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۱ 
g(j,1,s) = V(s);4-V(s)g,Bi “2,3 


where V(s)g, is the vector of components of V(s) correspond- 
i 


ing to the variables x, which are basic at B, (and recall 


that & 4 th 


4 İs the j 


column of A). Since V(s) = V! * s(V"-V') 
this expression can be written: 
ች » " -1 
¿(ASI ۴ (Vp. *s(V" p -Y'p1))B4 E 
: -1 ~1 / 
mesay < OBC .) + Be. Ne) BE an . : 
م‎ Nm 5 5 TU ታ፦" 3) (V 8, V و8‎ (23 5.4 | 
Substituting in the identity s = s* + (s-s*) give: 
P: " 08 =ጊ 1ه‎ 
Pis)” 72ء ۷۳)- 00009 سو سد‎ 5.1) 


-] 
ncs . - n <= 5 ۰ 
+(s-s*)[(V j ۷۲ (و‎ (V B 7" ,و‎ (5 54) 


— : : =1 
ኗ(3,3,5*)*(»-=*) ])۷۷ ج۷۷ ج۷۷)۔(۷۶۔‎ (3 "Ta; 


= e(j,1,s*)+(s-s*)2(j,i) 


Lemma V-9: g(j,i,s) and ¢(j,i) are defined as in the previous 
lemma. Let T be the set of all alternate bases B; solving 

for fixed s (in the interval O $ s $ 1) the problem: 

Max V(s)x 

subject to x in F 
Some basis B; in the set T has one of the following proper- 
ties: 

(1) e(3,1) É O for all j 


(2) s' » minimum -gl(j,i,s)/g(j,i) >0 
j/s(j,i)> 0 
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First it will be shown that there exists some basis in T 
which solves, for some positive e, the linear program: 

Max V(ste)x 

SuUbgeccs vo XIN F 
Select arbitrary small e and solve the problem above. Let T' 
be the set of bases solving this problem. If some member of 
T* is also a member of T, proceed to the second part of the 
proof. Otherwise, let Bi be some member of T! and notice 
that (from the reduced cost identity): 

g(j,i,s*e) = e(j,i,s)te g(j,i) for all j (*) 

Since B; solves the problem above, g(j,i,ste) = 0 for all : 
But since B; is not in T (hence does not solve Max V(s)x for 
x in F) g(j,i,s)> O for at least one j. Find a new value 
of e as follows: 


e! = 172 minimum (-e(3,1,8)/2(3,1)) 
j/eldsi,s) > 0 


new e = minimum (e', 1/2 of old e) 
e' is positive, because from (*), and the comment immedi- 
ately following, when g(j,i,s) is positive, g(j,i) must be 
negative. Then e is positive also, and is smaller than the 
previous e. Resolve the linear program for the new value of 
e, and repeat the above analysis. Continue to repeat until 
some basis in T is found to solve the probiem for some posi- 
tive e. Such a basis will ultimately appear, for the 


following reasons. First, the boundedness of F guarantees 


br = me 
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that each linear program has a solution. Secondly, e' is 
computed to be the smallest value of e such that the tableau 
associated with the associated basis satisfies the optimali- 
ty criterion. Since the new e is smaller than e', and every 
succeeding e is smaller yet, the basis giving rise to e' can 
never again appear ina later T;. But the number of bases 
is finite; eventually all bases not in T must be exhausted, 
and then the desired result must obtain. 
(second part) Let Bi be the basis found which solves both 
Max V(ste)x subject to x in F 
and 
Max V(s) subject to x in F 

The two optimality criteria this basis satisfies are: 

g(j,i,ste) " g(j,i,s) + e g( 3,1) $ O for all j 

e(j,i,s) £ O for all j 
Together these imply that one of the following conditions 
must hold: 

(1) g(3,i) £ 0 for all j 
(2) s' = minimum -g(j,i,s)/g(j,i) exists and is > 0. 
j/g(3,1,s) » 0 


QED 


Theorem V-10 ("Connectivity Theorem"): If B' and B" are two differ- 
ent bases whose tableaux satisfy the efficiency criterion, 


then there exists & sequence of neighboring bases, whose 





Proof: 


PII 


tableaux satisfy the efficiency criterion, linking tne two. 
Let A' and A" and x' and x" be the matrices of tight con- 
straints and the vertices corresponding to B' and B". By 
corollary V-3-2, there exist vectors k', k", L'>0, and L'>0 
such that k'A! = L'C and k"A" = L"C, and where when con- 
straint i is an inequality in the statement of F then 

E 11 = s 1 11 
kı' = O, or k = 0, as appropriate. For some such k', k”, 


L', L" define: 


0 Eike 
vo = T 
If V' = V" then B' and B" are alternate optimal bases to 


Max L"Cx subject to x in F 
By tne connectivity lemma (lemma V-7) the desired sequences 
exists. If V' f£ V" define V(s) = V! + s(V"-V') for scalar s. 
Notice that V(O) = V', V(1) » V", and 


V(s) = (sL" + (1-s)L')C. For s in the interval 0 e 


1, 

the vector multiplying C is strictly positive, so by theoren 

۷7-2 all solutions to "Max V(s)x subject to x in F" for fixed 

s in that interval are efficient. 

It is desired to find the bases which are optimal to the 

following linear program as s is increased from zero to one: 
Max V(s)x subject to x in F 

This is a well kmown single parameter parametric problem; the 


simplex may be employed to find the succession of optima 


which arise &s s is increased. 
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The bases so generated constitute the desired sequence. In the 
following, the method of construction of this sequence is explained, 
and it is demonstrated that only efficient bases are included, and 
that B" must ultimately be reached. Define ¢(j,1,s) and g( j,i) as 
in the statement of the reduced costs identity. 
Set i= 1, B' = By, and sy = O and commence the following: 

(Iteration beginning) 

75٦1613655۰1 Therefore, Pi ss) £0 for all Jo Recall‏ قرط 


that the reduced costs identity specifies: 
g(3,1,5:44) = g(j,i,s;) و‎ (54”=4+1)5(3›3) 


Several cases may arise: 
Case (a): One of the following two conditions obtains: 
(3) e(3,1) £ O for 211 j 


(2) S34] = minimum ل( )3 و‎ Si exists 
3/8(3,1)>0 


and is greater than or equal to one. (The 
determination of s,+,] in this fashion will be 
called the "subproblem"). 


In either of these events, B: 


; and B" are alternate optima to 


Max V(s.)x subject to x in F 
By the connectivity lemma (lemma V-7) the desired sequence exists. 
Case (b): The solution S34} to the subproblem above is > 54 
but less than l, Then apply the exist criterion of the simplex to 


any problem attaining the subdproblem minimum, and pivot as indicated 
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(the boundedness of F ensures that the exit criterion will be satis- 
fied by some basic variable). Label the new basis Bi+]» This basis 


is feasible and, since it solves 
Max V(si44)x subject to x in F 


for s3+ ] in the O-] interval, it is efficient. Of course, it neigh- 
bors B;. Set i = itl and commence another iteration. 
Case (c): The solution to the subproblem exists but = s;. 


In this event, there are alternate optimal bases to the problem 
Max V(s;)x subject to x in F 


By lemma V-9 at least one of these (call it Bz) has ons of the 


following properties: 


ዘቦ 


(1) se @ tan Oe Gor qe 15 


(2) Si+] = minimum (-z(j,t,s,)/glj,t)) exists and is > O. 
j/e(j.t)> 0 


By the connectivity lemma (lemma V-7) this basis can be reached by a 
finite sequence of neighboring efficient bases. Once reached, the 
situation will fall into one of the categories in case (a) or (b) 


and the remarks there apply. 


Notice that the sequence so constructed will terminate only when B" 
is reached. Moreover, the number of bases generated in the course 
of each iteration is finite. Also, s; is increased each iteration 


in such a manner that the basis used in its construction can not be 
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optimal for succeeding s Since the number of bases is finite, B" 


4° 
must ultimately appear. 


QED 


The following definitions and abbreviations will be employed in the 
statement of the algorithm for finding all efficient vertices to 
problem P and in Chapter VI. 
Problem set AP: This is the family of problems which deter- 
the efficiency of the bases neighboring a given basis. 
Let G be the reduced costs matrix for the basis whose 
neighbors it is desired to evaluate, and suppose that 
the dimensions of G are k by q, where k is the number 
of rows of C and q is the number of nonbasic variables, 
including sla variables, in ne constraints definings 
F. Let N be a k-vector, T a q-vector (with components 
T;), I the rank q identity matrix, and 1 a vector of 
ones. Problem set AP is the following linsar program 
for each i such that x; is nonbasic: 
Mind. 
set. GYM + IT = - Gt] 


M,T = 0 


LEB; List of efficient bases. (It is sufficient to keep 


only the indices of the basic columns). Bases are 





indexed "LEV;¡”. 


115 


An Algorithm for Finding ልጊጊ Vertices Efficient to Problem P. 


Step #1: 


Find an initial efficient basis; enter it on the 


LEB as LEB}. 565 3 = 13, 3 = 0. 


tep #2: 


Set j = jtl. If j > i, stop. Otherwise generate 


the tableau associated with LEB j ° 


Step #3: 
(2) 


(b) 


Do the following for each nonbasic variable xy: 
Ascertain which variables could leave the basis if 
the variable xy were to enter. The exit criteria 


ares 


d, = O and Yer 7 0, or 
9 دآ‎ (8 / sx) 
S/Ys,k > 0 


LEB is the list of the basic variables in the 
current basis. For each basic variable meeting 
the exis sriteria above, form a new list of basic 
variables by substituting x, for this exiting 
variable in the LEB 4. If at least one such new 
basis is not already listed on the LEB, solve the 
corresponding problem AP. If this problem has a 


solution, add the new unlisted bases to the LEB 


and increment i by l for each addition. 


When all nonbasic variables have been so considered, zo to 


step #2. 
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Theorem V-11: If x" is an efficient vertex to problem P, then some 


Proof: 


basis representing x" is on the LEB at the conclusion of the 
preceding algorithm. 

Let B' be the initial efficient basis found in step Fl. By 
theorem V-6, some basis representing efficient vertex x" has 
a tableau satisfying the efficiency criterion. Let one such 
basis be B". If B" f B' then by theorem V-10 there exists a 
sequence cf neighboring bases, whose tableaux satisfy the 
efficiency criterion, linking the two. The neighbor by 
neighbor search of steps #2 and #3 uncovers each of these 


bases, including B". 
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V-5 AN ALGORITHM FOR FINDING THE EFFICIENT MANIFOLD 

Not all] efficient solutions to problem P need be vertices; 
finding expressions for those which are not is the task of this 
section. 

Assume that all efficient vertices are known. It was 
remarked in section V-2 that facets may be expressed as convex com-~ 
binations of vertices of F; clearly, then, efficient facets may be 
expressed as convex combinations of efficient vertices of F. A set 
of efficient facets may be expressed in terms of groups of efficient 
vertices where it is understood that any convex combination of the 
vertices in each group is also efficient. This form of expression 
may be used to specify the efficient manifold provided that all 
efficient solutions are contained in efficient facets. The follow- 


ing lemma and theorem establish that proviso. 


Lemma V-12: If solution x is efficient to problem P then at least 
one constraint is tight at x. 

Proof: If problem P contains equality constraints, the conclusion 
is immediate. If not, suppose that x is efficient but no 
constraint is tight at x. From the statement of problem P 
there does not exist row vector L > O such that LC = 0; 
therefore, by Steimke's theorem of the alternative, there 
exists x' such that Cx' 2 O. Let K be the set of all j such 
that aj3x' 2 O (recall that the coefficients of the nonneza- 


tivity constraints &re included in A). 
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ማኝ minimum ( (b,-8,x) /a x") if K is not vacuous 
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bet x = + ex! 


Notice that Asx" = a.x + ۲و‎ 


There are two cases: 


£ : ፌ : ፡ : i 
(a) a¿x' = 0; then a,x" = b; so x" satisfies constraint i. 


(b) 8۳: 9; 


x" 


| 


as 


4 a;x +(minimum ( (b 5-2 5x) /a sx") ) ajx' 


17-7 ژ 


ዘኩ 


asx + ((54-84=)/ጾዩ4=')5ጎ4=፤ 


ዘኑ 


bi 
therefore constraint i is satisfied at x". 

Thenx is in Es 

But Cx" = Cx + eCx' = Cx because e > O, Cx! = O, and at 

least one component (Cx'); » 0. 


1 


Then is not efficient; a contradictions 


Theorem V-13: If point x is efficient to problem P then x is con- 
tained in an efficient facet. 

Proof: Let J be the set of indices of constraints tight at x (by 
the previcus lemma, J is not vacuous). Consider the facet 
ceneraved by une constraime set J; Ghis facet is efficient 
(by corollary V-3-2) if there exist vector L >O and multi- 


pliers k, (nonnegative when constraint i is an inequality) 





Jag 


such that 


LC = 7 eno 
i in 
But since x is efficient, such L and ks exist. 


QED 


Notice that every efficient vertex is some efficient facet. 


It is appropriate at this point to comment on some other 
conclusions about facets. It is a natural first assumption that 
facets containing efficient vertices must be efficient. This is not 
true. Nor is it true that a facet every extreme point of which is 
efficient is necessarily efficient (a counterexample is included in 
Appendix II). On the other hand, a facet containing an inefficient 
vertex can not be efficient. This useful result will be taken up in 
theorem V-15. A most valuable implication of this is that the defin- 
ing index set of every nonvacuous efficient facet must be contained 
in the defining index set of some efficient vertex. 

The reader will recall that "facet" was defined in associa- 
tion with a collection of tight constraints. This suggests the 
following useful concept: 

"Generic member": The "generic member" of a facet (or, 

occasionally, a “generic point" in a facet) will be 


defined to be any solution x satisfying as tight 
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constraints those and only those constraints which 
define the facet. 
The generic member nesd not actwally exist for all facets 

(that is, there may be no point in F actually satisfying those and 
only those constraints as equalities); nevertheless, whether it 
exists or not, it will be said to be efficient if (oy corollary 
V-3-2) there exist L 2 E > O and multipliers ki ( = 0 when constraint 
i is an inequality) such that for the index set J of defining con- 


straints 


d‏ 08د د1 


"Generically efficient": A facet will be said to be gener- 
ically efficient کل‎ 175 2enerie member 15 efficiıent, 
"Not generically efficient": A facet will be said to be not 
generically efficient if its generic member is not 
8۲ 17 ۰ 
In addition to generic members, facets may also contain 
other members (among them, perhaps, the vertices of the facets, 
which are also vertices of F) which satisfy additional constraints 
as tight constraints. Consider some such point satisfying the addi- 
tional tight constraints indexed by the set J'. Clearly if tre gen- 
eric member is efficient, then this solution is also; to show this, 
it suffices t edd k; = Ost نی وہ ۶۸و‎ Ug the LESE و وكا‎ iin J, 


satisfying the corollary V-3-2 condition for the generic point, to 
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which satisfies the efficiency conditions of that corollary for this 
other member as well. It follows that the efficiency of the generic 
member is sufficient to show the efficiency of the facet. 

Testing tne efficiency of the generic member can bs cast 
into a linear programming format by the by now familiar device of 
substituting M = L-1 into the condition cited; the condition beccmes 


Find M;,k i in J, such that 


12 


/; MC - ). ask; = - 10 


i in J i in Jd 


ir’ 


Ms 0 


1 


iy 


K: 


i; 7 O when constraint i is an inequality con- 


strand, 


If no feasible solution to this point exists, then no generic member 
of the facet is efficient, although nothing can be said about other 
members of this facet.  However--to preview theorem V-15--clearly 
there &re no efficient solutions to P whose tight constreints are a 
subset of J. 

in a problem with n + m variables and slack variables andm 
constraints, one can form zum _] different nonvacuous index sets. 


One way to find the efficient manifold is to find and test the 
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efficiency of eacn of the facets these index sets define. Having 
done so, the vacuous ones may be somehow weeded out, and the 
remainder represent the efficient manifold. The technique is simple, 
foolproof, and impractical for all but the smallest problems. 
Fortunately there are numerous shortcuts which can be employed. 

One, that only subsets of vertex defining index sets need be con- 
sidered, was mentioned earlier. The rollowing theorems motivate two 


more. 


Theorem V-14 (E - inference theorem): If J and J' are the con- 
straint index sets defining two facets, if f(J) is effi- 
cient, end if J' includes J, then f(J') is efficient also. 

Proof: The efficiency of the first facet implies there exist L= E>0 


andik: (nonnegative or inequality constraints) such that 


LC = 2. aK 
1 2 ال‎ 


Then 


۱۳۲ < Pa asks + 7 211 


1 1n ow 1 ل1‎ 
NO AnS 


Doms. = OS e in Ji! butonobiin J, sosthe) seneric member sof 


f(J') and hence f(J') itself are efficient. 


QED 
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Theorem V-15 (NGE-inference theorem): If J and J' are the defining 
index Sets of رت ا‎ and ٣٩ د ول‎ J includes J’, and if f(J) 
is not generically efficient, then f(J') is not generically 
efficient. 

ool: If Ghere does not exist L> 0 and k; ( = O when constraint i 


is an inequality) such that 


LC = 7 aik; + 2; asks 
i in J and J' ion dJ 
1 Aot and" 


then there can not exist L > ۵0 and k, (restricted as before) 


such that 


QED 


Fortunately, the E-inference theorem and the NGE-inference 
theorem can be used to considerably shorten the facet search pro- 
cedure outlined earlier. To begin with, the reader should be aware 
that it is not generally necessary to list all efficient facets in 
order to specify the efficient manifold. This is because often some 
efficient facets are contained in others (this is the case when the 
defining index set of the second contains that of the first). When- 
ever such cases occur, listing the large facet suffices to identify 


not only all efficient solutions contained in it, but also those 
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contained in the smaller facet as well. It follows that in the 
course of searching for efficient facets,when one is fond, its 
including facets (that is, those defining index sets containing its 
defining index set) can safely be ignored. Conversely, when a not 
genericaliy efficient facet is found, those facets including it (that 
is, whose defining index sets are included in the index set of the 
first facet) can be ignored. This is the principle alluded to 
earlier, wnen it was noted that only defining index sets contained 

in the defining index set of some efficient vertex need be con- 
Sidered. 

Unfortunately these twin principles work in opposite 
directions. When one is looking at facets of one constraint, the 
E-inference theorem is potentially quite valuable, because the 
opportunity for these facets to inelude others is high. On the 
other hand, the NGE inference theorem is worthless; no nonvecuous 
proper subsets of the list of tight constraints exist, so no effi- 
cient facets can contain these. If one's facet search progresses 
from large facets (with few defining indices) toward smaller con- 


. *. 


straints (with many defining indices) then the NGE-inference theorem 
remains useless, for all it can tell is what has already been pains- 
takingly discovered--that larger including facets are not generical- 
ly efficient. Conversely, a facet search which progresses from 
large index set facets toward small index set facets finds no value 


in the E-inference theorem. Starting with middlesized facets 


reduces the potential of both theorems. The upshot is that there 
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seems to be no a priori grounds for preferring any particular facet 
search above all others. The choice for a particular application 
must revolve about practical considerations including tne nature of 
the problem and the cambilities of the computer at hand. The effi- 
cient facet search procedure about to be presented is motivated by 
the anticipated virtue of a "self fathoming” feature (an opportunity 
to avoid useless ccmputations by structuring the sequence of opera- 
tions to take full advantage of the E-inference theorem) which it 
exhibits. 

The following definitions and abbreviations will be employed 


in the statement of the algorithm (and in the following chapter): 


"Problem Q" This is the problem which determines whether 
the generic member of a facet (hence the facet itself) 
is efficient. Let M be a k-vector, 1 a vector of ones, 
and k, scalar multipliers (nonnegative whenever con- 
straint i is an inequality constraint). "The problem © 
associated with set J" (or with facet f(J)) is to find 


M and the k; such that 


MC - 2. ask = -1C 
23٤1 


= Q when constraint i is an inequality con- 
straint. 
(This is a problem for phase one of the two phase 


simplex). 





LEM 


LEF.: 


LOF; 


An Efficient 


Step 
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List of efficient vertices. The information stored 
fom each vertex is the indices of whe constraints which 
are tight at that vertex. LEV; indicates the ith 
vertex stored on the LEV. 
The LEV is easily constructed from the complete 
LEB by ascertaining for each member of the LER 
which constraints are tight at the vertex repre- 
sented. The resulting vertex list is then checked 
for duplication (which will arise whenever several 
bases from a single degenerate vertex appear on 
the. ALBEDO 
List of efficient facets. The information stored for 
each facet is its defining index set. LEF; indicates 
the iP facet on the LEF. 
List of candidate facets. Facets are specified by 
their defining index sets; LCF; is the ith facet on 
this list. The function of this list will become clear 


once the algorithm is understood. 


Facet Finding Algorithm 

#1: Store the defining index sets of the efficient 
werbices. ineLEV, 41 tel, 2, cml. wet 5.- Bhe setoof 
indices of equality constraints. Attempt to solve the 
problem Q associated with 5. If a solution exists, 


place S on the LEF and stop. Otherwise set i = 0. 
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Step #2: Set i= itl. If i> U stop. Otherwise set J' = S; 
set J" = the defining index set of LEV, 
= (5, 1 ¿Loy - ipl Set LCF = the empty set. 

Step $3: Set i, = the index with the smallest subscript in 
J" - S which is greater than all index subscripts in 
J'- S. If there is no such index, go to step 76. 
Otherwisessetek es 

Step 74: Form the index set J' * i,. Determine whether 
this index set contains any index set on the LEF; if 
so, go to step #5. If not, attempt to solve the asso- 
ciated problem Q. If a solution exists, place this 
index set on the LEF; otherwise place it at the bottom 
of the LCF. 

Step $5: Set k= k+l. If kpr go to step #6; otherwise 
go to step #4, 

Step #6: If the LCF is empty, go to step #2. Otherwise 
remove the top index set from the LCF; call this set 


J'. Go to step #. 


This algorithm is clearly finite because the number of ver- 
tices is finite, the number of subsets of the defining index set of 
any verte is finite, andino sueset is generated more than once in 
the course of investigating the defining index set of a given 


vertex. 
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Example V-1: 

Suppose that LEV, of a given problem has the defining index 
set (S سن درد دسا‎ where iy indexes the jth constraint, and suppose 
that of the LEF for this problem, two facets Bin) and Dr )) 
have index sets contained in the defining index set of LEV). Then, 
substituting the values of the subscripts for the constraint indices 
--e.g. "3" signifies iz, and representing the problem Q associated 
with index set J by Q(J), the algorithm vrocesses the defining index 

sets arising from this vertex as follows: 

Q(S) is found not to have a solution. 

Q(S,1), 0(S,2), Q(S,3) are found not to have solutions; the 

setom S oT (2) (S5) are plaged on the ECKE. 

Q(S,4) nas à solution; (S,4) is placed on the LEF. 

k now exceeds r (which is 4) so the top member of the LCF, 

(S,1) is removed (and k is reset). 
Q(S,1,2) is found to have no solution; (S,1,2) goes to the 
bottom of the LCF. 

QUS.1,3) has à solution; (9,1,23) is placed on the LEF. 

Set (S,1,4) is generated but found to contain (8,4) from the 

LEF; 15 is discarded. 

(S,2) is removed from the LCF. 

Q(S,2,3) has no solution; (S,2,3) is placed at the bottom of 

the LCF. 

Set (S,2,4) is generated but found to contain (onu. 

(8,3) is removed from the LCF; (S,3,4) is generated but 


found to contain (S,4). 
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(S,1,2) is removed from the LCF; both (S,1,2,3) and (S,1,2, 
4) are found to contain members of the LEF. 

(S,2,3) is removed from the LCF; (S,2,3,4) contains a set on 
the LEF. 


The LCF is now empty; the algorithm passes to LEV»y. 


The relevant portion of the LEF has indeed been found; moreover, 


while there are 15 possible subsets of LEV,, only 8 attempted solu- 


tions of problems Q were required. 


Theorem V-16: Solution x is efficient to problem P if anà only if x 


Proof: 


is contained in some facet on the LEF constructed by the 
preceding algorithm. 

(If clause): A facet is placed on the LEF (in step #4) only 
if its corresponding problem Q has a solution; this is 
sufficient for tne efficiency of all solutions contained in 
the facet. (Only if clause): Let x be efficient to problem 
P. By theorem V-15 x is contained in some efficient facet. 
Let J" be the defining index set of one such facet and let J 
be the smallest subset of J” (there may be several--any one 
will do) containing S (the indices of the equality con- 
straints) and defining an efficient facet. If J = S then J 
is placed on the LEF in step #1. Otherwise, note that since 
f(J) contains x, it contains at least one vertex cn the LEV. 
Consider the first such vertex listed there; its defining 


index set contains J. In due course this vertex is taken up 
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for consideration in step #2. In step #4 thereafter a set 
containing S and the inequality constraint index in J with 
the lowest subscrivt is formed. Since this is the first 


vertex containing J, J (and of course any subset) can not 


ርኸ 


already be on the LEF. If this index set is J, then the 
corresponding problem Q is solved and J is placed on the 
LEF. Otherwise no solution to problem Q exists (guaranteed 
by the method of construction of J) so the index set is 
placed» on the LCF. In due courses it is recovered and the 
inequality index of next lowest subscript is added. If this 
is J, it is placed on the LEF; otherwise, it is placed on 
the LCF. The process continues av least unti" Jis CO 
structed in its entirety and placed on the LEF. 
QED 
The facet searcn procedure just outlined is reasonably easy 
to understand and straightforward to program, yet it does exhibit 
some computational niceties. For example, no problem Q associated 
with an efficient facet is solved twice; after its first appearance, 
the facet is listed on the LEF, which is checked before each problem 
Q is attempted. Moreover, no facets contained in other efficient 
facets are investigated; this is the "self fathominz" feature. On 
the other nand, storing the LCF may sometimes be a considerable 
burden, and &s written, nothing prevents repeated attempts to find 


solutions to problems Q corresponding to facets which are not 
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generically efficient. 

it @sppessible, by Wtroducing a™"slight™’conmplication, to 
remedy the last mentioned fault. In order to accomplish this, an 
additional operation must be inserted in step #4 between the LEF 
check and the attempt to solve problem Q. It is desirable at this 
point to check to see if the index set under consideration is con- 
tained in the defining index set of some previously investigated 
member of the LEV. If this is the case, then place the set on the 
LCF and go to step #5. In effect, the set has been inferred to 
represent a not generically efficiens facet, for if it ere effi- 
cient, it or a subset would have appeared on the LEF during con- 
sideration of a previous LEV. With this adaptation, the algorithm's 
handling of problems with many efficient vertices is expedited con- 
siderably. 

The goal of this chapter has been to develop a way to 
express the efficient manifold. It is now possible to do so. The 
efficient manifold is composed of one or more efficient facets; each 
of these may be expressed as the set of all convex combinations of 
its vertices. Vertices are easily associated with the facets con- 
taining then; a vertex is in a facet if the defining index set of 
the former contains that of the latter. In addition to being speci- 
fied in terms of defining index sets, each vertex has been identi- 
fied as a solution "x". This information was required for applica- 
tion of the exit criteria in the efficient vertex search, and for 


tne conversion of the LEB into the LEV. Defining the vertices of 
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eacn facet in terms of the soiutions x which they represent 
accomplishes a precise specification of the efficient manifold, the 


set of all efficient solutions to problem P. 
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CHAPTER VI 


ALTERNATIVE FACET FINDING TECHNIQUES 


11-1 GENERAL 

The previous chapter introduced a method for finding an 
expression for the efficient manifold. That method was comprised of 
two phases; an algorithm for finding the set of efficient vertices, 
and another for finding efficient facets. The "phase two" algorithm 
already presented is but one of many ways to go about finding these 
efficient facets. Tnree significant characteristics of the applica- 
tion of a given facet search procedure to a given problem are the 
number of problems Q attempted, the greatest lengths attained by the 
various lists kept in conjunction with the operation of the algor- 
ithm, and the complexity of the algorithm itself. These constitute 
measures of the running time of the algorithm, its storage require- 
ments, and its computer coding complexity. Though it is desirable 
to minimize all three, they are to some degree inherently conflict- 
inge The importance to be attached to any one of them will vary 
witn the circumstances of the user and the character of the problem 
at hand. There is no ideal facet search procedure apparent. 

This chapter discusses three alternative ways to carry out 
the second phase of the efficient manifold presentation method. The 


three approaches are: 
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(a) A facet search algorithm similar to phase two of the 


(c) 


method of Chapter V, except that the search progresses 
from facets of many constraints to facets of few con- 
straints. This algorithm (the "bottoms-up” algorithm) 
is presented in section VI-2. 

A facet search procedure combining the approaches of 
(a) abova and of the second phase of the Chapter Y 
approach. This algorithm (the "double-ended" algor- 
ithn) is presented in section VI-3. 

An approach to facet finding which requires less infor- 
mation to be stored in the coursewof the algorithm but 
which can be expected to require more subproblems 
(linear programs) to be solved than any other facet 
search procedure discussed. This approach (the "low 


storage algorithms") is presented in section VI-4. 
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VI-2 A “BOTTOMS-UP" EFFICIENT FACET FINDING ALGORITHM 

In the "bottoms-up" algorithm, the efficient facet search 
begins by considering facets with large numbers of tight constraints, 
and works toward those with few such. As in all of the facet search 
procedures, the search is founded upon the tight constraint sets of 
the efficient vertices--which, the reader will recall, are facets 
themselves. The following names and abbreviations will be used in 
the statement of this algorithm: 

LEV: The list of efficient vertices, each specified by the 
list of indices of the constraints which are satisfied 
as equalities at the vertex in question. 

LCF: <A list of candidate facets, specified by lists of their 
tight constraint indices. 

LNGF: A list of not generically efficient facets, similarly 
specified. 

LEF: The list of efficient facets, similarly specified. 

Q(J): The problem Q associated with the constraint index 
set J. 

Defining index set: A set of constraint indices. The 
defining index set corresponding to a vertex or facet 
is the set of indices of constraints satisfied as 
equalities. 

The members of each list will be designated by subscripting the sym- 
bol for the list--e.g., LEV, is the defining index set representing 


the yen vertex on the list LEV. 
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The function of the variable "flag" in the bottoms-up and 
double-ended algorithms is to keep track of whether or not a given 
defining index set under consideration has given rise to another 
index set, a subset of the first, which is stored on the list of 
candidate facets. If not, flag = 0 is the signal for the index set 


to be placed on the LEF. 


The bottoms-up algorithm: 

Step #1: Store the defining index sets of the efficient 
vertices in LEVi, i 7 1,2,...,U. Set S = the set of 
indices of the equality constraints. Attempt to solve 
the problem 9 associated with S. If a solution exists, 
viace S on the LEF and stop. Otherwise set i = 0. 

Step #2: Set i= itl. If i >U, stop. Otherwise set J" = 
the defining index set of LEV; 7 (8,1,,15,...,1.). Set 
J = J". Set LCF and LNGF = the empty set. 

Step 73: Set variable "flag" = 0. Find i,, the inequality 
constraint index with the lowest subscript in J which 
is higher than all indices in J” - J. If no such i, 

exists, go to step +6. Otherwise set k = u. 

Step #4: Form index set J' = J - 3... If J' is contained in 
8 9 ود لال‎ j < 1, or Centained in some LNGF 4, eo to 
step 75. If J' contains LEF; set flag = 1, place J' at 
the bottom of the LCF, and go to step #5. Otherwise 


attempt to solve problem Q(J'). If a solution exists, 
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set flag = l and place J' at the bottom of the LCF. 
Otherwise add J' to the LNGF (remove from the LNGF any 
sets contained in J'). 


Step #5: Set k= k + 1. If k Fr, go to step #4. 


Step #6: If flag = O and if J contains no set in the LEF, 
add J to the LEF and remove sets containing J. If the 
LCF is empty, go to step #2. Otherwise remove the top- 


most member of the LCF; call this set J; go to step oe 


Example VI-1: 
In this example the bottoms-up algorithm is applied to the 


problem of Example V-1. Recall that LEV) has the defining index set 
(S,i],i2,i3,i4) ana that the relevant portion of the LEF is 


የ(5,11,13) and £(S,ig)» Hereafter, subscripted indices will be 
represented by their subscripts alone; Q(J) indicates the problem Q 
associated with the set J; "Yes" will mean that this problem has a 


"in" indicates either that the 


solution; "No" indicates otherwise. 
set is a subset or superset of the following list; whether subset or 
superset will be clear from the context. "to" means the set is to 
be placed on the list indicated. When this is the LCF, placement at 
the bottom is understood. For example "(S,1,4) is in LEF" says that 
the set SGT ee, contains sane set on the LEF; "Q(S,4) Yes; to LCF" 
means that the problem Q associated with the set (S,i,) has been 


attempted and a solution is found, and that the set (S,i,) is to be 


placed on the LCF, at the bottom. Only the sets and problems Q 
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arising in step #4 of the algorithm will be listed conpletely; other 
information will be supplemented where clarity demands. The alzor- 
Jabmeproeeeeasc as follows: 

Q(S) No 

ND‏ ا 

MS, 20374) Tes: to ECF; "SE >= 1 

OCS. O) to LCF 

Be) Voss to LER 

0) co) 1562 77 

The LCF yields set J = (S,2,3,4); flag = 0 

SS LITE COO 77ء سن‎ 

OOO AES: TO LEP 

Q(S,2,3) No: to LNGF 

The LCF yields set J = (S,1,3,4)3 flag = 0 

Q(S,1,4) Yes: to LCF; flag = 1 

mo lo) Yes: to LOF 

The LCF yields set J = (S,1,2,4): flag = 0 

Q(S,1,2) No: to LNGF 

flag is still = 0 when k exceeds r: (6,1,2,4) to LEF 


LCF yields set J = (S,1,2,3): flag remains zero: to LEF 


EGF yields set J = (8,3,4): flag = O 
Q(S,4) Yes: to LCF: flag = 1 

(8,3) is in LNGF 

LCF yields J = (S,2,4): flag reset = 0. 


(S,2) in LNGF 
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flag remains = O when k exceeds r: (S,2,4) to the LEF: 
remove (S,1,2,4) from the LEF 

LCF yields J = (S,1,4): reset flag = 0 

(S,1) is in LNGF: flag = 0 and k exceeds r: (S,1,4) to LEF 


LCF yields J = (S,1,3): flag O 
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(S,1) is in LNGF: flag remains = 0: (S,1,3) to LEF (and 

remove (S,1,2,3)) 

LCF yields J = (S,4): reset flag = ۵ 

(S) is in LNGF: flag is still = 0: (S,4) to LEF (and 

remove (S,2,4)) 

The LCF is now empty: the algorithm proceeds to consider 
LEV, etc. Twelve different problems Q were attempted here (none 
more than once). 

$ will now be shown that the bottoms-up algorithm finds all 
efficient facets. 

Consider I, a set containing the equality constraint indices 
and an arbitrary subset of the inequality constraint indices of the 
defining index set of some LEV. Let J be the complement of I with 
respect to LEV,, where J * (j1,j2,...,j9) &nd jg2J444 for all n. 


Define the following sequence of index sets: 


ün s Eo JSt —.* jew ub Xndqk seb LEV) 


dioe uet سول‎ do * “an Je.1 
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This sequence will be called the “discovery sequence” for I from 


LEV. 


m; because it is the sequence of sets leading to set I in the 


course of investam@atime thewsubsets ofetnesdefinaaceindex=setwoof 
LEV. That is, suppose that I is generated in step #4. This set can 
only have been generated from the set Ją, by the following reason- 
ing. I is generated from some set by the elimination of one member; 
then this predecessor set must contain I and one other element--that 
is, one element of J. Suppose that element were something other 
than jj. Then ول‎ is not in this prior set, and the extra element 
which is, is less than .ول‎ But since only indices greater than or 
equal to i, may be eliminated, and since i, is greater than all 


indices of the LEV which are missing from the set in question (step 


#3), and since the extra element is less than the missing ài that 


Ca 


۰ ۰ ۰ ij - 
element can not be eliminated in step 74. By the same reasoning, J. 


must have been generated from J,_} 


1, and so on, until ول‎ the ver- 


tex's index set, is reached. 

Remark: The discovery sequence for set I is unique because there is 
only one increasing order of the set J. It follows that the 
bottoms-up algorithm generates each set 1 no more than once 
-~-thus it attempts to solve problem Q{I) at most once. 

Lemma Viel: If set I is placed on the LCF, then from that time for- 
ward I or a subset will be listed on the LCF or LEF or be 
under consideration as the set J in steps #3 to #6. 

Proof: When a set J on the LCF is removed, in step #6, then in the 


steps leading up to and including the next return to step 





Lemma 


Proofs 


#6, either some subset of J is placed on the LCF (flag = 1) 
or else (flag remaining = 0) J is placed on the LEF. Sets 

on the LEF are removed only when proper subsets of them are 
placed thergon. 


QED 


VI-2: Let LEV, be the first efficient vertex whose defining 
index set contains set I, which contains 5, the set of 
equality constraint indices. Problem Q(I) has a solution if 
and only if I or some subset of I is contained in the LEF 
when the bottoms-up algorithm concludes the consideration of 
the subsets of the defining index set of LEV. 

(If clause): Sets are placed on the LEF only after their 
problems Q are shown to have solutions. 

):6 مس‎ MÀ Assume that problem Q(I) has a solution. 
It will be shown that I or a subset of I must be on the LEF 
by the conclusion of consideration of LEV» 

Consider ihe sets J4,J5,..., in the discovery sequence for I 
fron LEV,. Since I is a subset of each, no prior LEV con- 
tains any of them. Also, by theorem V-14, solutions exist 
to each of their problems Q. In step +1, Q(S) is tested 
and if a solution exists, S 1s placed on the LEF: the 
desired result is obtained. Otherwise the algorithm per- 
forms the following operations (smong others) when LEV, is 


reached. 





Theorem 


Proofs: 
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Jo is generated and since it represents an efficient facet, 
it is placed on the LCF before step #6 is reached. Eventu- 
ally Jo is removed from the LCF, Jz is generated, and placed 
on the LCF, again before step #6 is reached. The algorithm 
proceeds to generate this sequence, some member of which is 
on the LCF whenever step #6 is reached, at least until the 
original set I is reached and considered. Of course, I, 
too, is placed on the LCF. By lemma VI-l, either I ora 
subset must be on the LEF at the conclusion of consideration 
of LEV... for the LCF is then empty. 


QED 


“ 


Vivo Solution x is erficient vo problem P if amdeonibyedf 

x is contained in same facet on the final LEF constructed by 
the bottoms-up algorithn. 

(If clause): Clearly x is efficient if it is contained ina 
facet on the LEF, for facets are listed thereon only if their 
prodlems Q have been shown to have solutions. 

(Only if clause): Assume that x is efficient. By theorem 
V-13, x is contained in some efficient facet f(J'). This 
facet is nonvacuous (it contains x) so J' is contained in 
the defining index set of at least one efficient vertex. 

Let LEV. be the first such on the LEV. By the previous 
lemma, J' or some subset is on the LEF at the conclusion of 


the algorithm. Let this set be J. Facet f(J) contains 
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سوت)5کسومسمفسووہ: ( نت 


The bottoms-up algorithm can be simplified considerably, at 
some cost to its efficiency, by eliminating all reference to the 
LNGF, in which case there may occur more than one attempt to solve 
problems Q for which solutions do not exist, and by simplifying the 
rules of formation of the sets J' in steps #3 and #4 to be simply 
“form all sets which may be formed by the omission of a single 
inequality constraint index in J", in which case some redundancy 
will occur on the LCF (in effect, the discovery sequence for a given 
set will no longer be unique). 

An embellishment which might pay for itself by reducing 
storage requirements would be to periodically purge from the LNGF 
any subsets of other sets on the LNGF. A convenient time to do this 


is when new sets are added to that list. 
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VI-3 A DOUBLE-ENDED FACET FINDING ALGORITHM 
The algorithm described here is a simple combination of the 
algorithms of the previous section and of section V-5. The facets 
containing each efficient vertex are investigated starting with both 
the largest (few defining indices) and the smallest (many defining 
indices) and concluding with facets of an intermediate dimension. 
The two halves of the algorithm can profit from one another by the 
mechanisms of theorems V-14 and V-15. This appealing feature is 
gained at the expense of considerable added complexity, including 
two lists of candidate facets. The two ares 
LCF1: The list of candidate facets used in the He portion 
of the algorithm, the part patterned after the algo- 
rithm of section V-5. Members of this list are not 
Poder ically wert volent. 
LCF2: The list of candidate facets used in the ¿na portion 
--the bottoms-up portion of the algorithm. These 
facets are generically efficient. 
With the excsption of the LCF, the other lists defined heretofore 
are required in this algorithm also. List members are indexad and 
specified as in the previous algorithms. A new function and two 
additional variables are also required. 
c(J) is defined to be the number of indices--the cardinality 
--of the set J. 
c' and c" are the cardinalities of the most recently 


investigated members of LCFl and LCF2. 
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The double-ended algorithm progresses through the possibili- 
ties at each efficient vertex by checking all candidate facets of 
cardinality one, then all candidates of the largest possible cardi- 
nality, then all facets of cardinality two, then all of cardinality 
one less than the largest, and so on. When ۵۱۲ + 1 = c” then all 
possibilities of any cardinality have been investigated; the algo- 


rithn moves on to the next efficient vertex. 


The Double-Ended Algorithm: 

Step #1: Attempt to solve the problem Q associated with the 
equality index set S; if a solution exists, store S on 
the LEF and stop. Otherwise store the defining index 
sets of the efficient vertices in LEV;, i = 1,2,...,U. 


Then set i = 0. 


Step #2: Set i = i + 1. If i> U, stope Otherwise set 
LCF1=LCF2=LNGF= the empty set. Set J" = the defining 
index set of LEV, = (e ió and let J = S, 

Jits Jand c ae) 

Step #5: Set c' = c(J). Find i, = the index in J" with the 
smallest subscript greater than all index subscripts in 
J-S. If i, does not exist go to step #6. Otherwise 
get - us 

Step #4: Form the set J' = J*iy. If J' contains some index 
set on ths LEF, go to step #5. If J! is contained in 


some prior LEV, or in some set on the LNCF, add J' to 
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the bottom of LCFl and go to step $5. Otherwise 
attempt to solve the associated problem Q. If a solu- 
tion exists, add*J' to the LEF; otherwise add it to the 
LCF1. 

Step #5: Set k= ktl. If عا‎ = r, go to step #4. 

Step #6: If LCFl is empty, go to step $11. Otherwise 
remove the top member of that list; call it J. If 
c(J) = c' go to step #3. Otherwise set c' = c(J). If 
c(J)*1 » c" go to step #11. 

Step #7: Set c" = c(J"'), flag = 0, and find i, = the index 
in J"' with the smallest subscript greater than all 
index subscripts in J" - J"'. If no such i, exists go 
to step #10. Otherwise set k= ue 

Step #8: Form the set K - J"'-iy. If set K is contained in 
some prior member of the LEV or in some set on the LNGF 
then go to step #9. If K contains some index set on 
the LEF, set flag = 1 and place K at the bottom of the 
LCF2. Otherwise attempt to solve Q(K). If a solution 
exists set flag = 1 and place K at the bottom of LCF2. 
Otherwise add K to the LNGF. 

Step #9: Set k= k+l. If k € r go to step 48. 


Step #10: If flag = O place” set K on the LEF. If LCF2 is 


“If a set to be placed on the LEF contains a set already 
thereon, the new set should be ignored. Sets containing a newly 
added set should be removed from the LEF. 
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empty go to step 711. Otherwise remove the top member 
of the 130۱20060۱۱ 129۵ ۲۰ If c(J"') - c" go to step * 
37. Obthermuse set c= c(J"'). If e” > ctrl goto 
step #3. 

Step 711: Place the remaining members of the LCF2 on the 


LEF. Go to step #2. 


۷1-2: Let LEV, = (S,1,2,3,4,5) and lst f(S,2) be the only 
efficient facet containing LEV. The notation of the previ- 
ous examples will be used. The double ended algorithm pro- 
ceeds as follows. 

Q(S) No 

Q(S,1) No: to LCF1 

05:2 Yes: to LEE 

6 Sy Nor’ toner! 

WS 4) No: to ECFi 

Q(S,5) Nos: to LCF1 

Now the bottoms-up portion of the algorithm commences: 
(S,2,3,4,5) is in LEF: to LCF2 

Q(S,1,3,4,5) No: to LNGF 

(SOIN S) 19/18 ۳۰ ها‎ 2 

(3 1998955) is in WES to LOF2 

(5,1,2,3,4) is in IF: to. LOEO 

Now the Chapter V portion, again: 


J= (StI) from CCFL. (8,1,2) Ts in LEF 





11 ር ር 1 
) ۸ dais ገር ሬን 10021 


CS تسا‎ sa ing CE to TEL 


J = (S,8) fran LCF1: (S,3,4) and (S,3,5) are in LNGF: both 
to LCF1 

J = (8,4) fram LCF1: ($,4,5) is in LNGF: to LCF1 

J = (S,5); no additional set is generated from this set. 


The algorithm returns to the bottoms-up portion of the algo- 
a Gi. 

gi! 05.2. 3.4.5). from LEF2. (S.34 5) is in 7 

(o Saa os Are all i تا‎ ۶ 

JS A Frauen. Beer 1د‎ (Ss. 172,2) are 
in LEF: 66 52 

T NM cU D رہ‎ UC. SO te ده‎ 

J'' 2 (S,1,2,5,4) from LCF2, but no subsets of this are cen- 
erated. 

J"! = (S,2,4,5) and it is found that c(J"') = c'; the algo- 
ritnm proceeds to step 11, finds that no additions to 
the LEF are required, and passes on to the next member 


of the LEV. 


Only seven problems Q were attempted in the course of this solution, 
thanxs to the early placement of an entry on the LEF and another on 
the LNGF. Bither the bottoms-up algorithm or the algorithm of 

Chapter V would attempt seventeen problems Q in the course of investi- 


gating this efficient facet. 
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This algorithm is clearly finite, for the following reasons. 
The number of vertices is finite, and the progression of the algo- 
rithm does not repeat vertices. The number of subsets of the defin- 
ing index set of any particular vertex is also finite. There are 
two enumeration procedures at work in the algorithm and neither 
repeats itself or the other; no subset of a particular vertex defin- 
ing index set is generated more than once in the process of investi- 
gating that vertex. Then the consideration of each vertex is finite; 
the number of vertex investigations is finite; the algorithm is 
finite. 

The discovery sequence discussed in conjunction with the 
bottons-up algorithm is pertinent to this algorithm also, in con- 
junction with its bottans-up portion, steps #6 - #9. Analagous to 
it is another discovery sequence operative during that portion of 
the algorithm patterned after Chapter V. This "downward discovery 
sequence" is constituted as follows. Let I be the set for which the 
sequence is to be constructed, and let it be composed of indices 
11,12,...,ip where i$ € ig. The enumeration process in the Chap- 


ter V portion of the algorithm generates the following sequence: 


Kj ij 
Ko = 17,19 


Kz = 31›32›13 
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Like the other sequence, once this downward discovery sequence is 
embarxed upon, there will always be one member of the sequence under 
consideration or on the list of candidate facets, in this case the 
LCF1, until some discovery implying the efficiency (in this case) or 
not generic efficiency (in the case of the downward discovery 
sequence) of the remainder of the sequence is made. The event which 
can terminate this sequence is the discovery of the efficiency of 
some facet in the sequence. It is immaterial whether the efficiency 
of the member is itself determined directly or whether the effi- 
ciency of the facet generated by some subset of a member of the 
sequence is established in conjunction with another portion of the 


algorithm. 


5۶560 5۱۹4۶ Solution x is efficient to problem P if and onlyeif x 
is contained in soma facet on the LEF at the conclusion of 
the double-ended algorithm. 

Proof: (If clause): If x is contained in a facet listed on the LEF 
then x is efficient, for facets are placed thereon only 
after their problems Q are shown to have solutions. 

(Only if clause): Let x be efficient to problem P. By 
theorem V-13, x is contained in some efficient facet. This 
facet is nonvacuous, since it contains x. Its defining 
index set J will be contained in the defining index set of 
some efficient vertex. Consider the situation at the end of 


consideration of the first such vertex on the LEV. 
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Case (a): c(J) £ ct. Consider the downward discovery 
sequence for J. Let J' be the first set thereon which 
defines a generically efficient facet. J contains J' so x 
is contained in f(J'). Since no member of the sequence 
prior to J' defines a generically efficient facet, each has 
generated its successor, which was placed on the LCFl. 
Ultimately J' was generated, found to describe an efficient 
facet, and placed on the LEF. 

Case (b): c(J) 3 c", Consider the discovery sequence 
for the set J. Every member of this sequence describes an 
efficient facet. Then every member has caused its successor 
to be placed on the LCF2. Then J can not fail to nave teen 
generated. Once placed on the LCF2 either it was placed on 
the LEF or a subset was placed thereon, either because no 
further efficient subset was generated, or because the set 
or subset was on the LCF2 when c'+] = c" obtained. 

Case (c): c' < c(J) < c" and the algorithm passed into 
step 411 when LCF1 was found to be empty in step $6. The 
remarks pertaining to case (a) apply; LCF1 can not become 
empty until some subset of J defining an efficient facet is 
found. 

Case (d): 6۱ > c(J)< c" and the algorithm passed into 
step #11 when LCF2 was found to be empty in step #10. The 
remarks of case (b) apoly; LCF2 can not become empty until J 


or a subset is found to be generically efficient. Then c(J) 
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must be = c" => a contradiction; case (d) can not arise. 


When a set has been placed on the LEF, it is subse- 
quently removed only if replaced oy a subset of itself. 
Phen at the conclusion of the algorithm some sst containing 
solutionsxəis «om the LEF. 


QED 


Notice that during the course of the double-ended algorithm, 
problem Q(J) for arbitrary J is attempted at most once. This may be 
seen from the following reasoning. Other than Q(S), problems Q are 
attempted only in steps #4 and #8. In each case the LEV is checked 
prior to the attempt to solve the problem, so a given problem can 
arise during the investigation of at most one vertex. Now consider 
the situation at the end of consideration of a vertex in which (ل)9‎ 
is investigated. Either c(J) $ c' or c(J) 2 c" but not both; then 
set J can not have arisen in both the Chapter V and bottoms-up por- 
ES oi the algorithm. The enimeratiom procedurs within each of 
these halves of the algorithm does not generate a given set more 
than once. 

This algorithm (and in fact any algorithm using the LNGF) 
Could begin with an initial List or not generically efficient index 
sets in the LNGF if the not generically efficient bases found in the 


course of the efficient vertex search are saved for this purpose. 
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VI-4 A LOW STORAGE EFFICIENT FACET FINDING PROCEDURE 

In the efficient facet search algorithms discussed hereto- 
fore, the lists LCF, LCF1l, LCF2, and LNGF were all used in various 
ways to avoid attempting to solve each given problem Q more than 
once» In addition, the LCF (or the pair LCF1, LCF2) guided the 
sequence of facets investigated, and indicated when to stop investi- 
gating facets arising in the context of a particular efficient ver- 
tex. The latter functions are the "self fathoming" effect mentioned 
briefly in Chapter V. The price paid for all of these computational 
conveniences is the need for storage of extensive lists. The facet 
search algorithms discussed in this section do not require the stor- 
age of lists other than the LEV and LEF. On the other hand, a signal 
advantage of these lists--the guarantee against multiple attempts at 
the solution of any particular problem Q--is lost. The one essential 
function which the omitted lists provided is the task of guiding the 
search through the possible alternative subsets of the defining 
index set of each efficient vertex in such a manner that no effi- 
cient facets are overlooked. In the algorithms about to be presented, 
this function is furnished by a sequencing rule which requires 
negligible storage for its operation. 

In the first algorithm here, sequencing is accomplished by 
associating each conceivable subset of the defining index set of the 
vertex in question with a nonnegative binary integer. Starting from 
zero, the integers are expressed in binary and this binary number is 


used to generate an index set. In this manner every possibility is 
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enumerated. (By contrast, the previous facet finders were able to 
avoid totally enumerating the index sets through capitalizing on the 
E-inference and NGE-inference theorems (V-14 and V-15) through the 
use of the LCF, or the LCFi and LCF2.) 

The binary number based sequencing rule operates as follows. 
To each index corresponding to an inequality constraint in the defin- 
joe index seu of the vertex under consideration 15 asgeciated a 
column in the binary number system. Columns are assigned starting 
from the right; thus if there are four indices of inequality con- 
streints in a given vertex's defining index set, then these would be 


11 
four's 


identified with the "one's column", the "two's column", the 
column" and the "eisht's column". The defining index set associ- 
ated with a given binary number contains two sets. The first is the 
Seu of indices of equelity constraintse Thessecond is» the set of 
those inequality constraint indices wnose columns in the binary num- 
ber contain a numeral 1, For exampie, if LEV, = UE igo doen) amd 3f 
ij is essociated with the one's column, ig with the two's column, 
and is with the four’s column, then the index set associated with 
the integer 6--in binary "110"--is (S,iz,ip). 

Bach time a new vertex is taken under consideration, the 
sets associated with the integers starting from zero are considered; 
this LEV is complete when the integers get too large to be expressed 
in binary with the columns to wnich indices have been assigned. In 


the example just given, 7 is the largest integer possible (and even 


this is not necessary, for it is the entire set) because 8 requires 
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four binary columns, and no index has been assigned to the fourth-- 


the eight's column. 


A low storage algorithm based on binary sequencing: 


Step 


Step 


Step 


Step 


#1: Store the defining index sets of the efficient ver- 
tices in LEVy,LEV>o, +... ,LEVy. Set i = 0. Attempt to 
solve the problem Q(S): if a solution exists, place 


the set (S) on the LEF and stop. Otherwise go to step 


i2: Set i = itl. If i exceeds U, stop. Otherwise set 
J" = the index set associated with LEF;. Associate 
(from the right) columns in the binary number system 
with inequality constraint indices in J". Set j = 1. 
#3: Ascertain whether the index set J associated with 
the number j expressed in binary contains an index set 
on the LEF. If not, attempt to solve tne associated 
problem Q. If a solution exists, add J to the LEF and 
remove from the LEF any index set containing Jd. 

#4; Set j = jtl. If j can not be expressed in binary 
with the columns to which indices are assigned, go to 


step $2. Otherwise go to step #3. 


Exemple VI-3: 


Consider the problem and notation establisned in example 


VI-1. (Recall that LEV, = (S,1,2,3,4); LEP includes f(S,1,3) and 


f(S,4)). The low storage algorithm with binary sequencing proceeds 
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as shovm in the following table: 





Number Binary Number: Set J 
(Baseyg) (44) (ig) (ig) (i4) 
1 1 (S,4) 
2 1 0 SS 
S 1 1 5 
4 1 0 0 ) 
5 1 0 (524 
6 1 1 0 (S 2.8 
7 1 1 Tas 20072) 
8 1 0 0 0 on 
9 1 0 0 1 (S,1,4) 
10 6. 0 (SAS) 
11 OI 
12 PO PEO (S, 130) 
13 1 1 0 1 O) 
14 je ጎዳ: በሙ 296922) 
15 1 1 1 رف اکا‎ 0 
16 


On the 
LEF? 


No 


No 


Yes 


Yes 


No 


Yes 


Yes 


Yes 


Bight different problems Q were attempted. 


Q(J)? 


No 


Yes 


No 


No 


No 


Yes 


No 
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Notes 


Q(S) 


to LEF 


to LEF 


too large to 


be repre- 
sented in 4 
binary 
columns; go 
to LEV, 
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inheorem VIS UU «x Ic eliicten, vorproulen P if and only if 
it is contained in some facet on the LEF generated by the 


algorithm just discussed. 


Corresponding to any possible facet is a binary number. 
This number must be generated in step #2. Proof of the theorem 
follows directly from this, by arguments now familiar from the 


proofs of each of the preceding algorithms. 


The following low storage algorithm employs a sequencing 
rule which provides the same sort of enumeration as is achieved by 
the algorithm of Chapter V. Since it omits the LCF of the earlier 
algorithm, it does not share that algorithm's efficient self fathom- 
ing feature. It does, however, have a crude self fathoming feature 
of its own based on the observation that if every conceivable facet 
whose defining index set has a given cardinality is efficient, then 


ever facet ol a larger cardinality must bB efficient also. 


A low storage facet finding algorithm with variable base sequencing: 
Step #1: Store the defining index sets of the efficient 
vertices in LEV,, i= 1,2,...,U. Set i= 0. Attempt 
to solve the problem Q associated with set S: ifa 
solution exists, place S on the LEF and stop. Other- 
wise set flag = QO. 
Step #2: Set i= itl. If i exceeds U, stop. Otherwise set 


J" = the defining index set of LEV.. Let n= the 
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mumber of indices in J"-S; associate each such index 
with an integer from 1 ton. Set j = 1. 

Step 73: Express j in the base ntl number system. Form the 
set J = S + every index whose associated numeral is in 
the base ntl expression for j. Ascertain whether ل‎ 
contains any set on the LEF; if not, attempt problem 
Q(J), and if a solution exists, place J on the LEF. 
Otherwise set flag = l. 

Step #4: Set j = the next larger integer expressible in the 
base ntl number system with the property that the 
numerals including zero in that expression are strictly 
decreasing from the right. If no such j exists, to 
step #2. If this j requires one more column in its 
base ntl exuression anan the ase and if flas ~ 0, 20 


to step #2. Otherwise set flag = O and go to step #3. 


Example VI-4: Let LEV, = (S,1,2,3,4) and let the efficient facets 
containing this vertex be £(S,1), £(S,3), and £(S,2,4). The algo- 
rithm just discussed would proceed with this problem as shown in the 
following table. Notation is as introduced in example Vi-1. "Skip" 
indicates that the base ntl representations of the integers con- 
cerned do nowt consist of numerals which are strictly increasing from 


the left. 
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Number 
(Base10) 


N? 


19-12 


Number 
(Bases) 


ነ 


he 


13 


23 


24 


Set J 


584 


On the 


LEF? 


Yes 


Yes 


Yes 


Yes 


No 


Yes 


Does Q(J) have 


& solution, if 
attempted? 


No 


Yes 


No 


Yes 
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Remarks 


Q(S) 


to LEF 


TOLLE 


skip 


skip 


to LEF 
skip 


All acceptable 
two digit base. 
4 2 
numoers have 
been shown to 
have efficient 
sets J; pro- 
ceed to LEV, 


Ten sets have been enumerated and six problems @ have been attempted. 
Were the previous low storage algorithm to be applied, it would also 
attempt six problems Q, but would enumerate all sixteen possible sets. 
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The proof of the completeness of this alzorithm is analagous 


to the proof of the completeness of the first low storage algorithm. 
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VI-5 COMMENTS AND SUMMARY 

The author has no computational experience with the facet 
search algorithms of this chapter. Nevertheless, it is perhaps well 
to summarize the expected attributes of these approaches. 

* The efficient facet search procedures of Chapter V and of 
section VI-2 are both straight forward in operation and 
reasonably simple to program for the computer. The 
choice between the two would depend most heavily on 
whether large (few tight constraints) or small (many 
tight constraints) facets are expected. 

* The low storage procedures are the simplest algorithms of 
the lot, both conceptually, and in programming terms. 
In terms of the numbers of linear programs to be 
solved, they can generally be assumed to be inferior to 
the other approaches. This is because only the most 
rudimentary precaution against multiple solutions of a 
given problem Q is used. In addition, the number of 
comparisons required may occasionally be quite large, 
because every conceivable alternative subset may be 
generated and inspected. 

* The double-ended facet finder has the greatest complexity 
and the greatest variety in its protections arainst 
duplication and solution of unnecessary lirear ٢ 7 
Wille it requires morê lists alld more conbuter 07 7٦ 


than the other approaches, its total storage requirements 
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might compare favorably with those of the Chapter V or 
the bottoms-up algorithms. Vihether the number of prob- 
lems Q attemoted favors this algorithm, the bottoms-up 
algorithm, or the algorithm of Chapter V depends on the 
character of the problem at hand. Example problems 


favoring each have been found. 


15 summary, the author is not at this time able to demon- 
strate the clear superiority of any of the efficient manifold 
searches cited in this chapter and the last. Apparently there are 
countless other possible facet search procedures. For exemple, an 
interesting variation on those of this chapter would be a bottoms-up 
low storage approach. Moreover, any of these procedures can be modi- 
fied in numerous ways (e.g., any one which employs the LNGF could 
choose to truncate it when it reaches a certain length). It would 
appear that the characteristics of all of these approaches are 
sufficiently diverse that any one of them could be at a disadvantage 
on sane occasion. In choosing one for implementation (or in design- 
ing a new one), the analyst should be guided by the expected magni- 
tude and complexity of the intended application, by the computer 
available, and perhaps even by the character of the linear program- 


ming packages available for solving the attendant subproblens. 






: ایس لہا۔ — 


Y መጪ LT و‎ 


1 5 اما 
ern | -— 2‏ 7 
0 ما SS‏ — 
E‏ یت تہ سس | غا 
SE a 4..ብ] ۱‏ 
جب ጋጋ‏ 
a‏ 
on ieee e ጩጨ ወ 8‏ = 


ር: © 
کک کاس‎ be Y -Í © 





















ا .د سا 
= — 
ےہ - . | 
ፕ 6‏ — ——- 


CHAPTER VII 


THE SCHOOL OF HEALTH SERVICES EQUILIBRIUM PLANNING MODEL: RESULTS 


VII-1 A REVIEW OF TEE EXPRESSION OF THE EFFICIENT MANIFOLD 

Before discussing the efficient solutions to the problem of 
Chapter IV, the manner in which the set of efficient solutions is to 
Bemexipecscame.) be prici iy revid or tae Senelit of the reader 
who has skipped Chapters V and VI. 

Imagine that it is desired to express all of the pairs 
(xj1,x2) which, if plotted, would fall in the cross Hatched region of 


77 76 IV. 
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One way to accomplish this is to list the corners ("vertices") of 
the shaded area, remarking that these and any points which can be 
expressed &s weighted sums of these, with nonnegative weights adding 
to 1, are in the shaded region. Such weighted sums are called "con- 
vex combinations". For example, in Figure V, the pair (3,2.5), 
which is labeled A, is a weighting of the vertices (2,3) and (4,2); 


yn fact, 


(SARS). = REN CEES Ca) 


Similarly the point B, representing the pair (11/3,3) is 


(ASS) =D TS) + 1/3 لش‎ 
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The set of efficient solutions, that is the “efficient mani- 
fold", to the School of Health Services equilibrium planning model 
is expressed in this manner--by specifying its vertices. These par- 
ticular results illustrate an occasion when the vertices must be 
organized into several sets. The totality of all convex combina- 
tions of the members of each such set is called a "facet". Members 
of different sets should not be used when forming convex combina- 
tions, for otherwise the result may not be efficient; i.e., may not 
fall within a facet. For example, Figure VI shows two (hatched) 
efficient facets together with the vertex sets generating them. 
Vertices B and D belong to both sets; that is, they are each con- 
tained in both facets. Incidentally, a consequence of the theory of 
Chepter V is that whenever an efficient manifold is constituted of 
more than one facet, then there will be enough vertices common to 
two or more facets that the efficient manifold is connected, as this 
one is. Vertex A belongs only to the first facet and vertex B 
belongs only to the second facet. Point E is a convex combination 
of A and C; notice that it is not a member of either efficient facet 


and is therefore not efficient. 
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VII-2 THE EFFICIENT MANIFOLD 

When the efficient manifold finding technique of Chapter V 
is applied to the problem expressed in section IV-4, 108 efficient 
vertices are found. These are arrayed in two efficient facets. 
Facet #1, containing 88-efficientevertices, is the set of all effi- 
cient solutions having additional clinical supervisor availability 
as their common limiting constraint. Facet #2 , containing 52 effi- 
cient vertices, is the set of all efficient solutions having faculty 
office and dry laboratory availability as their canmon limiting con- 
straints. It happens that 32 vertices are common to both facets. 

One way to characterize the vertices is by the values which 
the problem variables assume at them. This is what was done when 
the 108 count just mentioned was found. The model contains 17 vari- 
ables including slack variables (a slack variable is a quantity 
representing the amount of underutilization of the resource on the 
right hand side of one of the “less than or equal to” constraints). 
۱۲۰ 15 8190 0۵۵55019 to characterize the vertices according to the 
criteria values achieved at them; this is the naturai way for the 
decision maker to view them. There are only five criteria, and 
these are determined from only four variables. In particular, the 
slack veriables and the variables Y, through Y,, wnich are concerned 
with office and classroom utilization, do not figure directly in the 
annual numbers of graduates from each curriculum or in the annual 
surplus of the school. Consequently efficient vertices which differ 


only in their values of these variables are, from the standpoint of 
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the criteria identical» This phenonenon,.is so widespread in these 
particular results that among the 108 different efficient vertices 
only fifteen different criteria scores patterns are represented, 
Nine of these patterns belong to Facet #1, including four which also 
belong to Facet #2. Including the latter four, there are ten 
patterns in the second facet. Table I is a listing of these 
patterns according to whether they belong only to the first facet 
(Group A), to both facets (Group B), or only to the second facet 
(Group C). The values of the first four criteria have been rounded 
off to the nearest integer. These are, of course, the annual num- 
bers of graduates from the programs indicated. Since it has been 
ሠ... that there will be mo dropouts, and since each curriculum 
| 5 38 Wo year duration, these figures may also be considered to 
represent the annual matriculations or half of the annual enroll- 


1 


menes in the curricula represenved. 

The reader is reminded that the surpluses shown in Table I, 
Mee measurements accordlns to the fifth criterion, are variable 
instructional surplus. That is, only those costs and incomes which 
are Judged to De proportional to enrollments or faculty size are 
1 11:64. in order to find the total surplus corresponding to any 
of these it would be necessary to add to the figure snown another 
firure representing the sum of all other incomes such as grants 
unreiated to the enrollments of the School, minus the sum of all 


other expenses, such as the salaries of the faculty specialists. 


But notice that the relative impact of two different errollment and 
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Group À: 

A 

B 
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D 
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F 

G 

H 

I 
Group C: 
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Vertices in facet #1 alone: 
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M emt o 
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Vertices in facet #1 and facet #2: 
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47 
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Vertices in facet #2 alone: 
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Surplus ($000) 


-821 


-449 
-579 
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-575 
-920 
-884 
-657 
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Minimum Maximum 


Annual graduates: 


Total 69 159 

cdesioa b 0 69‏ و 
program‏ 

H. Associate 0 69 
and Nurse 
Practi- 
tioner 
programs 

H. Services O آ0‎ 
Mgmt. 
program 

Environmental 0 TOZ 
Hygiene 
program 

Instructional $-1,128 30 


Surplus ($000) 


TABLE LI 
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staffing programs can bs obtained simply by comparing the surpluses 
shown in the table. For example, pattern I, with a variable surplus 
of $-1,111,000 is $662,000 more expensive per year than pattern FL, 
with a variable surplus of 2-449,000, regardless of the magnitude of 
the non-proportional incomes and expenses. 

It happens that the maxima and minima of the values attained 
beyethe criteria ewer the set of all efficient. solutions are achieved 
at vertices. The maxima and minima for the criteria of this problem 
and the maximun and minimum values of the sum of the annual gradua- 
tions, or matriculations, are shown in Table II. These figures 
indicate the approximate dimension of the efficient manifold. 

Any convex combination of any of the patterns in groups A 
and B in Table I corresponds to an efficient solution in Facet #1. 
Similarly any sen combin=tionof some or all of thewpaUuberns in 
groups B and C is an efficient solution in Facet #2. Together these 
constitute the set of all efficient possibilities. Care should be 
taken to avoid mixing patterns from the first and third groups; E 
combina tens are not necessarily efficient. 

The following examples illustrate the use of Table I to gen- 
erate additional efficient matriculation patterns. Consider first a 
straightforward weighted averaging of members of Facet #1. Imagine 
that a decision maker has reviewed groups A and B-and decide&- to 
investigate the average of patterns © and H. He calculates vas 


follows: 





IRE 0, 03) 7,064) =01/2 of pattern C 


۱ ۳ ۰ ۲ 0 coy 0,-8lze)- 472 ofmpattern H 


il 


(34,31,43, 4, -581) 


That is, an equal averaging of these patterns yields an efficient 
solution in which 34 health assistants, 31 health associates or 
nurse practitioners, 45 health services managers, and 4 environ- 
mental hygienists are matriculated and graduated annually, and which 
has an annual variable instructional budget cost of $581,000.  Per- 
haps on examining tnis solution, the decision maker dscides that he 
would prefer more environmental hygienists; he might add a dash of 


pattern G, as follows: 


87 61 43133, 4, 2581 


il 


5 OF the previous result 


ti 


+ PROA OO S9) .2 of pattern G 


= (oS 34 291 nom) 


and so on. 

One may also experiment graphically with these possibilities. 
1 "۱5۸26 ۰۳ و‎ the criteria levels of patterns J amd P from group ۵ 
in Table 1 are marxed on ths left and right axes respectively. 
Notice that there are two scales shown; one for the numbers of 


raduates of the various programs, and the other for tne annual 
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variable surplus. Corresponding criteria for uhee two patterns Mavs 


boen joined by straicht lines and unese are laceled to sud whieh 
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criteria they represent. The intersections of these links with any 
vertical line drawn between the two vertical axes corresnond to the 
scores attained by some combination of the solutions J and P. For 
example, reading from the top down, the dashed line represents a 
pattern having 70 health services managers, 60 environmental hygien- 
ists, a surplus of roughly $-580,000, 30 health assistants, and no 
health associates or nurse practitioners graduating annually. Like 
patterns J and P, this new pattern is a member of Facet #2; there- 
fore it is efficient. 

If the decision maker wishes, he may transfer the results of 
one such graphical combination to a new axis, complement it with yet 
another pattern, and combine the two. A series of such successive 
combanetioms is» illustrated in Figure VIII. 

In addition to these camposite matriculations patterns, 
there is other interesting information provided by the efficient 
manifold. A perusal of the values of the problem variables at the 
efficient vertices reveals that no case exhausts the supply of 
demonstration laboratories or bench laboratories which was postu- 
186360 in Chapter IV. In fact, the lowest surplus capacities 
reported for these two laboratory types were 8 hours per week of 
demonstration laboratory availability, and 13 hours per week of 
bench laboratory availability. The reader will recall that no reli- 
able estimates of the true availability of bench or demonstration 
laboratory space was avaiiable, but that a reasonable minimum of ons 


full tims laboratory of each type was assumed. These resulis show 
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that if only this bare minimum were to be provided, the demonstra- 
tion laooratory would stand empty at least one day weekly, and the 
bench laboratory would be empty, at a minimum, almost two days 
weekly. 

The variable Y,, which represents the amount of usage of 
medium classrooms for undersized classes, 1s present at every effi- 
cient vertex; however, itis minimum value is only 4; that is, small 
classes are held in medium classrooms at least 4 hours weekly. When 
a small class occupies a large room 1% does so quite cemfortably, 
regardless of how much too large the room is. The reverse is not 
true. In allocating available classroom space to large, medium, and 
small classrooms, the greater evil is to err in the direction of 
rooms which ars too small. This posivive Y¿ throughout the efficient 
manifold confirms that this error nas been avoided. 

It is perhaps interesting to nove the computational charac- 
teristics of this solution. The author's computer code, which is an 
ilexpsrtly vrosmızmed example of une alsoritaajof Chapter Y, 
required 8 and 1/2 minutes for execution on an IBM 7094 computer. 
1702 problems "Q" were attempted, and the list of candidate facets 


attained a maximum length of 74 facets. 
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CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 


VIII-1 THE EFFICIENT MANIFOLD FINDING PROCEDURE 
The problem addressed in this thesis is to find the set of 


efficient solutions to 


(Maximize) F(x) 


subject to G(x) £0 


where F(x) and G(x) are linear vector valued functions of the vector 
x, a member of Euclidean n-space. A straightforward, finite pro- 
cedure (or donne SO has been revealed, and several warıaulonsson the 
basic method have been presented. The method and its several 
versions all rely in principle upon ordinary linear programming 
notions and in execution upon Simplex algorithm procedures. The 
efficient manifold finding procedure has been illustrated by appli- 
cation to a problem of planning in higher education, and the charac- 
teristics expected of the alternate facet finding procedures have 
been noted. 

The successful application of the efficient manifold finding 


procedure of Chapter V to the model of Chapter IV demonstrates that 
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the algorithm does indeed perform as anticipated. In Chapter VII it 
was established by example that the efficient manifold so presented 
can indeed provide helpful information. This success notwithstand- 
ing, certain improvements in the approach may be anticipated. 

Tne reader may have noted that the efficient manifold of 
Chapter VII is particularly simple. Mucn more involved results will 
occasionally arise from problems nearly as small, and quite moderate 
problems may have efficient manifolds which are beyond the capability 
of the decision maker to handle when presented in the artless 
manner illustrated. Clearly, better styles of presentation and 
better ways to cope with involved efficient manifolds are in order. 

Quite a number of variations on the basic efficient manifold 
finding procedure were suggested in Chapter VI. Exploration of these 
and other possibilities should continue; considerable streamlining 
may be possible. Two intriguing possibilities not previously men- 
tioned” are the followings” First, perhaps some sort of accelerator 
might be developed for use at degenerate vertices. The number of 
bases associated with a degenerate vertex grows monstrously as the 
number of excess zeros increases; depending on the circumstances, 
many or all of these may meet the efficiency criteria. The algo- 
rithns presented must investigate them all, just 8s if they repre- 
sented distinct vertices; perhaps this could be avoided. Secondly, 
amore complete merger of the vertex finding and facet finding 
phases of the approach might prove advantageous. Such an approach 


might involve making all calculations from a single tableau, a sort 
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of extended "problem Q" tableau (see Chapter V) in which every con- 
straint gradient is present but of which certain such gradients are 
ineligible for introduction into the basis in any given problem. It 
will be recalled that a vertex is a facet of cardinality m--which 
suggests a way to find an efficient, though not necessarily feasible, 
vertex directly frem this tableau. 

One of the particularly interesting and useful capabilities 
or goal programming is its ability to optimize several objectives 
seriaily--e.g., a second criterion might be optimized subject to the 
condition that the optimality of the first criterion is maintained. 
Such an option would be a welcome addition to the techniouss of the 
efficient manifold presentation method. The idea is to be able to 
eliminate a variety of multiple optimality which may arise when some 
variables are not represented in the natural criteria. Such a case 
arose in the model of Chapter IV; a mere 15 configurations of the 
1 51118685 Olzsınverest were represented by 108 distinct efficient 
vertices. It 1s possible to introduce an essentially arbitrary cri- 
terion in which the extra variables are featured, but if this new 
criterion is accorded equal status with the others then in addition 
to eliminating some redundancy it may also introduce new vertices 
which are efficient in the sense that some or all of the original 
criteria may be improved at the expense only of this secondary 
measure. Clearly these new solutions are quite as objectionable as 
the redundant ones which were eliminated. A preemptive priority 


procedure along the lines of the one used in goal prograrmine might 
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perhaps be used to produce the desired effect without introducing 
these unwanted vertices. 

Finally, there may be other good ways to go about finding 
efficient manifolds. It is possible, for example, that some benefit 
may be had from the knowledge of which vertices neighbor others. 


This information is immediately available from the lists of variables 


basic abeecachpivenevertexm 
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VIII-2 APPLICATIONS TO PLANNING IN HIGHER EDUCATION 

The last chapter reported on the application of the method 
of this thesis to a model of the School of Health Services. The 
Dean of that School found this small study to be rather enlightening; 
durinz the formulation of the model, certain previously unappreci- 
ated aspects of the School's operations were deservedly highlighted. 
Chief among these was the significance of the supoly and expense of 
part-time clinical supervisors. The nature of the efficient mani- 
fold which was revealed seems to confirm the importance of the 
availability of these particular faculty. In addition to providing 
early identification of the presently existing limitations to the 
school's growth, the efficient manifold is expected to be of value 
to tne staff of the School in preparing funding requests. 

In recent years the demand for analytic tools for institu- 
tional planning has markedly increased. This trend is a response to 
many pressures, including a sharpened appetite on the part of the 
state and federal governments for quantitative documentation of 
needs and performance. The chief impediment has been a vexing 
inability to satisfactorily treat with the outputs of higher educa- 
tion. A great deal of attention has been paid to this area--one 
bibliography lists almost two hundred references*--but the problems 
remain. A major difficulty seems to be that the products of educa- 
tion enjoy an ample inventory of menuinely incommensurable aspects. 
Impending developments and recent work in multicriterion problem 


solving may yet enable institutional planners to cope with this 
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troublesome diversity. The author believes that the efficient mani- 
fold presentation method holds particular promise because of its 
generous accommodation to the difficult circumstances of the 


decision maker. 
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FOOTNOTES 


acl. John H. and Robert D. Lamson, An Annotated Bibli- 
ography of Literature Relating to the Costs and Benefits of Graduate 
bamea tion, ane Council Of Graduate Schools, Washington, D.C., 1972. 





APPENDIX I 


In section V-3, following theorem V-6, it was observed that 
not all bases for a degenerate efficient vertex need have the 
property that for their respective reduced costs matrices G there 
exists vector L> O such that LG £ 0. The following example illus- 


trates such a case. 


Consider the following problem: 





(Max) E 0 ^ x]. 
| 
| | 
| و0‎ ۴ X3. 
3 = = E aL . , Slack Variables 
Subject to| 1 1 1/2: | 24 | 2 2 
Bea T 
| | | | 
0 1 0 | | ×3 ፡ 76 
| 35 لے‎ 
> 
X1,X2,Xa3 = 0 
The solutaon x' = )11 0 رر‎ is a degenerate vertex. Consider 


the basis B comprised of the columns corresponding to x,,x4,x5, and 


x73; from this, the reduced costs matrix G is developed, as follows: 
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Clearly there is no L>O such that LG = 0. Nevertheless this vertex 
is efficient, as will now be demonstrated. The constraints exactly 
satisfied by xl are the first, second, and third, and the nonnega- 
tivity constraint for xz, the relevant form of which is -xg £ 0 


The matrix A of the gradients of these constraints is 


I — 2 
1 0 O 
A = 
1 80 0 
| 
lo o -i 


For vectors k = (1,0,0,0), which is nonnegative, and L = (1,1,1/2), 
which is strictly positive, the condition kA = LC of theorem V~3 is 


met; tne point is efficient. 
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APPENDIX II 


In section V-5 it was remarked that a facet every extreme 
point of which is efficient ismot itself mecessarily efficient. 


The following example illustrates this point. Consider the problem 


seme ee 

(Max) | -10 l 1 | xy 
| 

l -10 1 و2‎ 

1 1 -10| 5۴3 


Subject toll E 1] x} | 3 (Slack variable: Xa) 


a a Ce  ፡ 
e 
CA 


X15X2,Xz = 0 


This problem is graphically represented in the following figure. The 


feaszible regions Gre etrenedson in the positive ortnant. 
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Consider the facet f(x,), the set of all feasible solutions for 
which the first constraint is tight. In the figure, this is the 
triangular face toward the veiwer. This facet has three vertices; 


they are 


EN o | o] 

| | ۱ 

o!‏ | 1 و 

| | 

xi = | | x? = | | x? 7 8 
0 ! |1.1 

| | | So 

| 

0 0 LO: 


Each of these will be shown to bs efficient. At xi, the first con- 
straint and the nonnegativity constraints for x5 and xz are tight; 


then the matrix of coefficients of tight constraints A is 


3 1 1 
A = lo S ٢٢ 
0 or ¿El 


Let k = (1,0,99) and L = (1,1,10); these are nonnegative and posi- 


tive respectively. Then 
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RAS GOS) | اھ وھ‎ oo = (152998) 
| 0 oe e E 
B | =10 1 1 
= (1,1,10) | ۳ -10 لژو‎ +7 
1 


and, by theorem V-3, x is efficient. Similarly, az x^, 
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and letting k = (1,0,99) and L = (1,1,10) we have kA = LC. 
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Lectin 7 21 99) and L= (MALO خم تور‎ = LC, 
Nevertheless, f(x ) is not efficient, for there are no k 3 0, L > O 


such that 
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Thus, for example, the feasible solution x* = (1/3,1/3,1/3,0), which 
is in f(x4) and satisfies no constraint other than that defining 


this facet, is not efficient. 
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